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1.  INTRODUCTION  & SUMMARY 


The  possibility  of  attaining  high  optical  extraction  efficiencies  from  lasers 
operating  on  the  bound -free  transitions  of  rare  gas  halide  excimers  formed 
in  h’ ’h  pressure  electrically  excited  mixtures  is  currently  attracting  much 
interest^^^"^^^^.  So  far,  the  most  successful  exploitation  of  this  scheme  has 
been  the  development  of  the  KrF  laser^^^  (22)^  extensive  literature^^^^ 
on  the  kinetic  processes  and  their  rate  constants  is  presently  available.  The 
present  program  has  been  directed  toward  experimental  and  theoretical 
investigation  of  the  fundamental  physical  mechanisms  of  the  KrF  laser,  with 
the  objective  of  demonstrating  high  volumetric  extraction  (~  50  J/liter)  at 
high  efficiency  (~  107c). 


In  order  to  demonstrate  the  potential  of  excimer  laser  systems  for  high 
optical  extraction  at  high  efficiencies,  a large  volume  laser  device  was 
constructed.  The  operating  characteristics  of  this  device,  which  shall  be 
described  in  detail  in  Section  2 , were  chosen  in  such  a way  as  to  provide 
maximum  flexibility  for  investigating,  in  addition  to  the  KrF  system,  other 
promising  molecules.  It  is  shown  how  the  basic  geometry  developed  from 
various  constraints  imposed  by  the  requirements  of  electron  beam  and 
discharge  excitation  of  the  low  impedance  medium.  Under  nominal  conditions, 
the  electron  gun  operated  at  a beam  voltage  of  3Z0  KeV  for  current  densities 
of  10-20  A/cm^  and  1000-2000  cm^  apertures.  Pulse  lengths,  typically  one 
microsecond,  required  the  development  of  discharge  circuits  capable  of 
driving  significant  energy  through  the  inherent  circuit  inductance  over  these 
times.  In  this  regard,  both  fast  capacitive  discharge  and  pulse  charged 
discharges,  designed  to  circumvent  a portion  of  the  circuit  inductance,  were 
developed.  The  resulting  circuitry  was  highly  successful  providing  for  a 
repeatable  and  reliable  source  at  full  rated  energy  and  charge  voltages. 

(l)-(30)  references  appear  at  the  end  of  Section  1.  (pages  2-4). 
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In  order  to  optimize  the  performance  of  the  KrF  laser  and  to  understand 
its  scaling  to  high  power  and  high  volume,  it  is  necessary  to  understand  the 
kinetic  and  pumping  processes  which  are  responsible  for  the  foimation  and 
quenching  of  the  excited  species.  In  Section  3 , the  details  of  a comprehensive 

theoretical  model^  ^ ^ will  be  desc’*ibed.  This  model,  which  consists  of 

an  extensive  computer  analysis  of  the  coupled  set  of  equations  for  molecular 
kinetics,  electron  kinetics,  radiative  extraction,  and  external  driving  circuit, 
was  constructed  from  a more  general  computer  code  which  automatically 
synthesizes  the  complete  coupled  analysis  by  translating  symbolic  reactions 
into  computer-coded  equations.  The  general  synthesis  code,  which  i?  applicable 
to  a broad  class  of  electrically  excited  laser  systems,  makes  i".  possible  to 
obtain  a complete  coupled  analysis  of  the  most  complicated  kinetic  scheme 
containing  an  arbitrary  number  of  reactions  and  interacting  species  with  virtually 
no  effort.  The  code  synthesized  for  analysis  of  the  KrF  system,  which  involved 
a reaction  scheme  of  ~ 80  kinetic  processes  has  been  successful  in  experimental 
comparisons  over  a v/ide  range  of  operating  parameters  for  the  scale-up 
and  other  experimental  devices.  Under  optiinum  conditions,  the  code  predicts  a 
specific  optical  output  of  ~ 35  J/liter  with  a 6.5%  extraction  efficiency. 

Section  4 contains  the  experimental  results  and  a comparison  with  the 

theoretical  analysis  for  a typical  set  of  conditions.  Due  to  nonuniform 

electron  beam  energy  deposition  in  the  direction  of  propagation,  both  total 

extracted  energy  and  discharge  enhancement  were  severely  limited.  As  a 

result,  experiments  were  conducted  with  3 and  10  liter  extraction  volumes  in. 

order  to  optimize  specific  energy  extraction  and  total  output  energy,  respectively. 

Under  optimum  conditions,  total  energy  output  was  90  J for  a 10  liter  volume, 

and  for  the  3 liter  volume,  specific  extraction  was  22  J/liter  for  beam  current 
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densities  on  the  order  of  20  A/cm  . Efficiency  in  both  cases  was  about  10%. 


^ ^ S.  K.  Searles  and  G.  A.  Hart,  Appl.  Phys.  Lett,  27  , 243  (1975);  Appl. 
Phys.  Lett.  79  (1974), 
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2.  DEVICE  CHARACTERISTI  S 

Electron  beam  and  discharge  excitation  of  large  volume,  gaseous  laser 
media  can  be  accomplished  only  through  a limited  set  of  geometries.  In  the 
present  case  of  high  pressure,  rare  gas -halogen  mixtures,  it  is  advantageous 
to  have  the  electron  beam  colinear  with  the  discharge  current  in  order  to 
increase  the  beam  range  which  would  otherwise  be  severely  limited.  The 
additional  requirement  of  a moderately  long  gain  medium  dictates  that  the 
optical  axis  be  aligned  at  right  angles  to  the  electron  beam  direction.  These 
constraints  ha\  o led  in  the  present  instance  to  a device  which  utilizes  a large 
aperture,  20  x 100  cm^,  e -beam  propagating  across  about  10  cm  of  gas  to 
a discharge  anode  of  similar  cross  sectional  area.  Various  other  system 
design  parameters  were  chosen  to  accomodate  the  use  of  high  pressure  and 
corrosive  gases  and  the  need  for  a low  inductance  discharge  circuit. 


The  * .>'?.l  laser  system  is  shown  in  a conceptual  rendering  in  Figure  1. 

The  high  pressure  laser  plenum  is  in  the  form  of  a half-cylinder  situated 

between  two  box  like  structures.  The  structure  to  the  left  of  the  plenum 

houses  the  electron  beam  power  supply  with  a 50  kJ  storage  capacity,  capable 

2 2 

of  sustaining  360  kV  beams  at  current  densities  of  20  A/cm  over  a 1000  cm 
aperture  for  0.  5 |i  s.  The  structure  to  the  right  contains  the  discharge  power 
supply  capable  of  100  kV  pulses  for  1 p.  s into  the  0.25  ohm  plasma  discharge. 
Other  peripheral  equipment  includes  a gas  handling  system,  seen  at  the  lower 
edge  of  the  drawing,  and  a shielded  electronics  diagnostics  room  at  the  upper 
edge.  These  components  are  described  in  detail  in  the  following  paragraphs. 


2.  1 Electron  Gun  . In  order  to  generate  the  large  area,  relatively  high 
electron  beam  current  densities  required  of  theKrF  system,  the  cold-cathode, 
field  emission  technology  was  adopted.  Such  devices  make  use  of  the  electrons 
produced  in  dense  plasmas  formed  around  sharp  cathode  edges  under  the 
action  of  intense  electric  fields.  The  electrons,  accelerated  in  the  generating 
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electric  field,  pass  through  a thin  metal  foil  which  serves  as  the  anode 

plane.  The  parent  plasma,  expanding  from  its  place  of  origin  at  the  cathode, 

approaches  the  anode  at  velocities  on  the  order  of  5 x 10  cm/sec  and  upon 

contact  terminates  the  pulse.  Since  the  background  pressure  is  on  the  order 

of  lO'^  mm  Hg,  the  residual  gas  plays  no  role  in  these  processes.  It  has 

been  shown  that  the  diode  behavior  is  space  charge  limited  and  therefore 

(31) 

follows  the  Child -Langmuir  law  , 


J ~ 


di 

d^(t) 


where  J is  current  density,  V is  electrode  voltage  and  d is  the  anode-cathode 
separation.  The  major  tradeoff  in  cold  cathode  devices  is  cherefore  between 
current  density  and  pulse  duration.  For  the  present  device  the  parametric 
range  for  beam  voltages  up  to  360  kV  is 


J = 22  A/cm^,  400  ns  pulse  length 

2 

12  A /cm  , 800  ns  pulse  length 
8 A/cm^,  1.2  |Jis  pulse  length 

2 

The  desired  beam  area  was  20  x 100  cm  , however,  depending  upon  current 
density,  magnetic  pinching  effects  in  the  diode  vacuum  section  greatly  limited 
the  attainable  emission  area.  A 3-stage  Marx  generator  was  used  for  high 
voltage,  high  current  power  generation  since  flexibility  was  required  in  the 
choice  of  operating  parameters.  It  was  further  required  that  the  current 
rise  time  be  on  the  order  of  100  nsec  and  this  was  achieved  with  a capacitive 
peaking  circuit  in  parallel  with  che  Marx  generator  output. 

The  electron  gun  is  therefore  composed  of  three  basic  elements:  (1)  the 
high  voltage  power  supply  in  the  form  of  a dc  charged,  three  stage  Marx  bank; 

(31) 

R.  H.  Huddlestone  and  S.  L.  Leonard,  "Plasma  Diagnostic  Techniques", 
Academic  Press,  1965,  page  119. 
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(2)  the  emitting  cathode  and  its  associated  corona  ring  and;  (3)  the  high 
pressure  foil  support.  These  items  are  indicated  pictorially  in  Figure  2 and 
will  be  discussed  below  in  detail  followed  by  a presentation  of  typical  perform- 

u . . (32) 

ance  characteristics 

Electrical  Design,  "he  electrical  system  consisted  basically  of  the  Marx 

bank  and  its  attendant  high  voltage  trigger  generators  and  crowbar  switch;  the 

output  peaking  capacitors  and  series  peaking  switch;  and  the  delay  trigger 

generators  and  control  circuitry  located  in  a console  external  to  the  device. 

All  of  the  high  voltage  components  were  immersed  in  insulating  oil  within 

a completely  sealed  metal  enclosure  to  prevent  excessive  rf  radiation.  Energy 

was  coupled  from  the  Marx  power  supply  to  the  gun  cathode  through  three 

feedthroughs  in  order  to  minimize  nonuniforniities  in  the  current  supplied  to 

the  cathode  surface.  Of  critical  importance  was  the  design  of  the  current 

feedthroughs  in  vacuum,  since  high  voltage  breakdown  was  a major  concern 

at  the  intended  gun  voltages.  As  a result,  a comprehensive  analysis  of  the 

(32) 

geometry  dependent  electrical  characteristics  was  undertaken  . Other 
elements  of  the  electrical  design  were  not  as  sent  ve  since  insulating  oil 
and  simply  excessive  stand-off  distance  could  be  relied  upon.  Metallic  walls 
surrounding  the  cathode  corona  ring,  for  example,  were  maintained  at  a 
distance  at  least  twice  that  from  the  ring  to  the  foil  ground  plane,  and  this 
arrangement  was  highly  satisfactory. 

The  crucial  elements  of  the  feedthrough  design  are  shown  in  Figure  3 
where  the  contours  indicate  equipotential  points  obtained  from  a computer 
calculation  of  Laplace's  equation.  The  areas  of  greatest  concern  are  along 
the  insulator  surface  and  at  the  insul'  ,or-conductor-vacuum  triple  point. 
Calculation  of  the  electric  field  along  '".he  insulator  indicates  a maximum  field 


(32) 


The  electron  gun  was  provided  by  Systems,  Science  and  Software  of 
Hayward,  California  who  performed  the  critical  design  study  as  well. 
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Feed  through 
Insulator  — 


Figure  Z Electron  Gun  Pictorial  Schematic 


Figure  3.  Feedthrough  Equipotential  Profiles 


intensity  of  57  kV/cm  for  400  kV  applied  to  the  cathode,  end  since  the 
predicted  breakdown  field  intensity  is  196  kV/cm,a  safety  factor  of  3.4  is 
realised.  The  triple  point  field  is  also  comfortably  low  based  on  design 
experience. 

The  gun  circuit  elements  were  specified  after  analysis  of  the  equivalent 

circuit  shown  in  Figure  4.  The  Marx  circuit  inductance  was  estimated  at 

about  1.  0 [iH,  therefore  it  was  expected  that  the  beam  voltage  risetime  would 

exceed  the  100  ns  limitation,  and  indeed  this  was  the  case.  In  order  to 

achieve  a shorter  risetime,  the  peaking  capacitor,  C^,  and  series  switch,  S^, 

were  added.  In  this  latter  case,  switch  was  set  to  actuate  after  reached 

the  desired  peak  voltage.  The  circuit  was  analyzed  for  various  values  of  C 

2 ^ 

and  an  initial  charge  voltage  of  400  kV.  A current  density  of  10  A /cm 
(20  kA  total)  was  assumed  and  the  nonlinear  dependence  of  current  density 
on  voltage  and  the  collapsing  anode-cathode  spacing  were  taken  into  account. 
The  results  are  shown  in  Figure  5 and  indicate  that  the  desired  risetimes 
are  achieveable  with  slight  overshoot  in  the  voltage,  dependent  on  the  choice 
of  C^.  These  data  were  reproduced  experimentally  in  their  essential  details 
as  outlined  below  in  the  discussion  of  gun  operation. 

With  regard  to  the  actual  circuit  components,  each  stage  of  the  Marx 
bank  was  composed  of  four  capacitors  in  a series /parallel  configuration  for 
a per  stage  capacitance  and  voltage  capability  of  1.85  pF  and  133  kV, 
respectively.  Therefore,  total  capacitance  and  maximum  erected  voltage 
for  the  three  identical  stages  was  0.  617  pF  at  400  kV,  at  a total  energy 
storage  of  50  kJ.  The  peaking  capacitor,  C^,  was  composed  of  2 banks  and 
each  capacitor  bank  was  formed  from  five  capacitors  in  series  for  a total 
bank  capacitance  of  40  nF  and  a stand-off  potential  of  500  kV. 

A mechanical  crowbar  was  included  at  the  output  of  the  Marx  bank  which 
disengaged  at  the  point  of  full  charge.  This  was  found  necessary  to  prevent 
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Figure  4.  Electron  Gun  Equivalent  Circuit 
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Cathode 


application  of  voltage  to  the  cathode  under  conditions  of  false  triggering 
which  otherwise  would  cause  erractic  erection  of  the  Marx  voltage  below  its 
intended  level  leading  to  arcing  in  the  vacuum  section  and  possible  foil 
rupture. 

The  gun  voltage  was  measured  with  a capacitive  divider  attached  directly 
at  the  output  of  the  Marx  bank  and  provided  a consistent  and  accurate  monitor 
of  the  gun  condition.  An  attempt  was  made  to  monitor  gun  current  at  each  of 
the  feedthroughs  and  at  the  Marx  ground  connection.  In  the  former  case, 
three  independent  Pearson  current  probes  were  installed,  one  surrounding 
each  of  the  three  current  feedthroughs.  However,  due  to  the  necessity  of  an 
extreme  high  voltage  stand-off  distance  between  the  feed  connection  and  the 
grounded  probe  case,  extraneous  signals  were  introduced  and  that  method  of 
measurement  was  abandoned.  The  ground  connection  monitor  was  found  to 
provide  an  accurate  measurement  of  total  current  but  was  not  useful  in 
diagnosing  various  gun  malfunctions  and  was  also  abandoned.  The  gun  voltage 
monitor,  correlated  with  Faraday  cup  measurements,  provided  the  only 
reliable  gun  diagnostic  tool. 

The  console  housed  all  the  controls  necessary  for  operation  of  the  gun. 
This  included  gauges  and  setting  valves  for  the  switch  SF^  pressures,  the 
delay  trigger  generator  for  actuation  of  the  start  and  stop  circuitry,  vacuum 
gauges  to  determine  the  condition  of  the  diode  section,  the  Marx  bank  charging 
voltage  indicator  and  various  interlock  mechanisms. 

Mechanical  Design.  Both  the  cathode  and  pressure  window  support  were 
central  to  the  gun  design  since  they  directly  determined  the  available  current 
densities  and  spatial  uniformity. 

_4 

The  cathode  was  composed  of  7.  6 x 10  cm  thick  tantalum  ribbon,  spot- 
welded  to  a stainless  steel  sheet  base  fixture  which  was  then  attached  to  the 
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face  of  the  corona  ring.  The  ribbons  projected  about  1 cm  above  the  support 

sheet  and  adjacent  ribbons,  running  the  1 M length  of  the  plate,  were  spaced 

approximately  1 cm  apart.  The  entire  emitting  surface,  constructed  in  this 

2 

manner,  was  25  x 100  cm  in  area.  Two  basic  types  of  cathode  were  employed 

one  providing  for  high  current  densities  on  the  order  of  20  A/cm^  and  the  other 

, 2 

providing  for  current  densities  up  to  about  12  A /cm  . Curvature  away  from 
the  anode  at  the  center  of  the  high  current  cathode  was  deemed  necessary  in 
order  to  increase  the  current  density  at  the  cathode  edges  and  so  minimize 
the  counter-effect  of  beam  magnetic  pinching.  By  the  Child- Langmuir  law 
presented  above,  the  central  current  density  would  necessarily  drop  below 
that  at  the  edges  as  a result  of  increased  anode-cathode  separation.  This 
was  found  to  be  the  case  and  will  be  discussed  below. 

The  diode  corona  ring  and  specially  j’  aped  field  grading  rings  at  the 
base  of  the  plastic  feedthroughs  were  made  of  polished  aluminum  and 
swabbed  with  diffusion  pump  oil  in  actual  operation.  Apparently  the  thin 
film  of  oil  greatly  reduced  the  tendency  to  arc  and  therefore  increased 
general  machine  reliability. 

The  cathode  structure  and  corona  ring  were  mounted  on  moveable, 
vacuum  tight  shafts  which  projected  through  the  plastic  feedthroughs  into 
the  Marx  tank.  This  allowed  variation  of  the  anode-cathode  distance  through  a 
mechanism  external  to  the  vacuum  vessel  without  destroying  the  diode  vacuum 
integrity. 

The  pressure  window  was  composed  of  a complex  sandwich  of  screens 

and  foils  as  shown  in  Figure  6.  In  their  path  from  the  cathode  surface  to  the 

-3 

high  pressure  cell,  the  electrons  are  seen  to  traverse  a 1.8  x 10  cm 

thick  aluminum  ground  foil,  a vacuum  drift  space  between  the  support  webs, 

-3 

a support  screen,  a 5 x 10  cm  thick  pressure  bearing  foil  and  finally  the 
arc  protection  screen. 
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Figure  6.  Foil  Support  Structure 
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The  ground  foil  was  necessary  to  provide  a uniform  ground  plane  for  the 
cathode  since  in  its  absence  arcing  would  occur  to  the  sharp  support  webs. 

This  effect  was  aggravated  by  the  requirement  of  high  current  densities  which 
necessitated  a small  anode-cathode  separation.  The  ground  foil  also  served 
to  absorb  low  energy  electrons  generated  during  the  voltage  rise,  thereby 
preventing  unnecessary  heating  of  the  second  foil  under  stress  by  high  pressure 
gas. 

The  support  web  structure  was  developed  to  provide  the  'greatest  open 

area  consistent  with  operating  pressures  up  to  5 atm  and  the  expected  beam 

incidence  angles  due  to  magnetic  deflection.  As  a result,  the  design  called  for 

deep,  widely  spaced  webs  with  a superimposed  support  screen.  The  support 

webs  were  formed  in  a ladder  configuration,  machined  from  a solid  plate  of 

aluminum.  All  the  webs  were  0.  8 cm  thick  and  spaced  so  as  to  provide  for 

2 

open  squares  of  49  -m  area.  Those  webs  ruruiing  longitudinally  down  the 
1 M extent  of  the  ladder  were  3 cm  deep  2.nd  those  in  the  2C  cm  transverse 
direction  were  2 cm  deep. 

The  support  screen  was  stamped  from  . 16  cm  thick  stainless  steel  in  a 
pattern  which  allov/ed  for  0.6  cm^  open  squares  and  8 x 10  ^ cm  wide  minor 
ribs.  The  resulting  configuration  of  support  webs  and  support  screen 
resulted  in  a geometrical  beam  transmission  factor  of  65%,  Of  course,  the 
actual  beam  transmission  factors  were  reduced  somewhat  from  this  value  due 
to  the  angular  trajectories  of  electrons  striking  the  support  ladder.  This 
effect  is  taken  up  below  in  consideration  of  beam  uniformity. 

.3 

The  5 X 10  cm  thick  high  pres  sure- vacuum  interface  foil  was  formed 
from  aluminum  alloy  #5056-H19,  since  this  afforded  the  highest  yield 
strength  for  the  foil  temperature  rise  expected.  The  maximum  pressure 
called  for  in  the  experimental  program  was  4 atm  absolute  and  under  this  condition 
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the  foil  and  its  support  structures  performed  satisfactorily  under  electron 
beam  operation. 

Finally,  a screen  was  placed  over  the  foil  but  separated  from  it  by  1 cm, 
?nd  served  the  function  of  providing  a uniform  ground  plane  opposite  the 
discharge  anode.  Secondarily,  arcs  generated  in  the  discharge  process  were 
terminated  on  the  screen  thereby  preventing  foil  rupture.  The  screen  was 
70%  transmissive  and  was  woven  from  0.  l6  cm  diameter  stainless  steel  wire. 

Diagnostics.  The  most  important  aspects  of  gun  performance  were 
current  density  and  deposition  uniformity.  In  order  to  determine  these 
properties  as  accurately  as  possible,  several  diagnostic  techniques  were 
employed.  With  regard  to  current  measurement,  a 0.  635  cm  diameter 
Faraday  cup  was  used  to  measure  temporal  behavior  and  a check  on 
absolute  current  was  made  with  a Pearson  current  probe  and  found  to 
correlate  within  10%^^^^  Overall  spatial  uniformity  was  measured  with 
electron  dosimetry  film  (cinemoid)  and  calibrated  with  the  Faraday  cup 
which  provided  excellent  spatial  resolution  owing  to  its  small  diameter. 

The  Faraday  cup  consisted  of  the  innermost,  0.635  cm  diameter  aluminum 
collector  surrounded  by  a pyrex  insulator  and  an  outer  stainless  steel  sleeve, 
connected  to  the  signal  cable  shield.  Crucial  to  proper  operation  was  the 
application  of  insulating  tape  over  the  probe  end  which  prevented  charge 
buildup  on  the  collector  from  shorting  either  to  the  shield  or  to  the  foil 
through  the  plasma  generated  in  air  by  the  electron  beam.  The  signal  was 
transmitted  directly  to  the  proper  attenuators  and  a 50  ohm  termination  at  the 
oscilloscope.  With  this  technique,  spurious  signals  generated  by  ElvU  were 
kept  to  approximately  an  oscilloscope  trace  width. 
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A check  on  the  Faraday  cup  accuracy  was  performed  with  the  apparatus 
shown  in  Figure  7.  The  entire  output  aperture  of  the  beam,  just  beyond  the 
foil  was  blanked  off  with  the  exception  of  a 2 cm  diameter  hole  in  the  center 
which  allowed  access  to  a small  portion  of  the  beam.  A plastic  insulatoi 
and  aluminum  collector  were  fitted  into  the  hole  and  the  collector  was  returned 
to  ground  through  the  Pearson  current  probe.  The  Faraday  cup  was  fitted 
to  the  same  blank  off  plate  to  provide  simultaneous  measurement.  As 
mentioned  above,  the  absolute  current  correlated  within  lOfo,  with  the 
Faraday  cup  reading  slightly  low,  and  no  difference  in  the  two  temporal 
waveforms  was  discernible. 

Spatial  plots  of  the  beam  output  aperture  were  made  with  electron 
dosimetry  film  and  simultaneously  with  the  Faraday  cup  over  a limited 
spatial  extent.  It  was  found  that  the  spatial  gradient  determined  from  the 
film  with  a transmission  densitometer  was  within  25%  of  that  measured  with 
the  Faraday  cup. 

Electron  Gun  Performance.  Typical  oscilloscope  traces  are  shown  in  Figure 
8.  Figure  8a  indicates  Marx  voltage  and  beam  current  measured  with  the 
Faraday  cup  at  the  foil  pressure  window,  in  the  absence  of  the  capacitive 
peaking  circuit.  The  voltage  is  shown  to  rise  to  roughly  half  its  maximum 
value  before  the  output  current  starts,  the  delay  being  due  to  the  requirement 
for  substantial  cathode  voltage  for  emission.  Total  voltage  risetime  is  200  ns, 
whereas  the  current  risetime  is  about  300  us.  Beyond  the  300  ns  point  in  the 
current  pulse,  the  slope  is  still  positive,  however,  this  further  increase  is 
attributable  to  the  decreasing  effective  anode-cathode  distance  as  the 
emitting  plasma  propagates  from  the  cathode  blades  to  the  anode  plane.  The 
current  waveform,  as  indicated  by  the  small  diameter  Faraday  cup,  shows 
only  slight  modulation  superimposed  on  the  basic  sinusoidal  shape.  Three 
Faraday  cups  separated  by  7 cm  in  the  output  window  major  direction  and 
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Figure  7.  Current  Density  Calibration  Apparatus 
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Figure  8.  Typical  Gun  Performance  Traces 


21 


operated  simultaneously,  showed  the  same  basic  temporal  shape  but  no 
correlation  in  the  modulation.  Therefore,  the  slight  noise  was  not  due  to 
common  mode  signals  and  was  spatially  randomized,  at  least  over  distances 
of  7 cm. 

Figure  8b  shows  both  Marx  voltage  and  resultant  output  current  with 
the  peaking  circuit  connected.  In  this  case,  the  voltage  risetime  is  about 
100  ns  with  a 20%  overshoot  beyond  the  average  value.  This  is  reflected 
in  the  fast  rising  current  pulse  and  leads  to  more  severe  modulation  which 
is  further  accentuated  by  the  plasma  propagation  mechanism  indicated  above. 
A series  of  Faraday  probes,  operated  as  outlined  above,  showed  the  same 
relative  oehaviors. 

Overall  uniformity  measurements  were  made  with  electron  dosimetry  film 
which  was  shown  to  be  sensitive  to  electron  current  but  insensitive  to  x-rays. 
Bleaching  of  the  film  was  therefore  related  directly  to  the  electron  current 
density  and  measurement  of  this  effect  with  an  optical  transmission 
densitometer  afforded  a direct  indication  of  beam  current  gradients.  Typical 
profiles  obtained  in  this  manner  with  the  film  in  contact  with  the  foil  are 
shown  in  Figures  9 and  10.  Figure  9 indicates  the  beam  current  distribution 
along  the  central  major  axis  of  the  output  window  for  the  high  density  cathode. 
The  accentuated  emission  at  each  end  is  due  to  the  slight  cathode  curvature 
as  mentioned  above  and  enhances  beam  uniformity  deposition  in  the  high 
pressure  gas  at  these  areas  which  would  otherwise  show  severe  end  gradients 
due  to  gas  scattering.  Similar  results  were  obtained  with  the  low  density 
cathode  at  central  current  densities  of  10  A/cm^.  In  this  latter  case,  the 
increase  in  end  intensity  was  not  as  pronounced. 

Figure  10  indicates  the  transverse  beam  uniformity,  scanned  from  the 
center  ci  the  foil  to  the  top  across  the  output  aperture  minor  axis. 
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Characteristic  of  the  high  density  cathode  I is  a severe  gradient  in  this 

dimension,  enhanced  by  the  presence  of  the  deep  supporting  webs.  The 

2 

fall-off  in  current  from  20  A /cm  at  the  center  to  roughly  half  this  value  at 
the  left  edge  of  the  web  is  likely  due  to  the  interception  of  those  electrons 
approaching  the  right  side  of  the  web  at  large  entrance  angles  due  to  vacuum 
magnetic  field  pinching.  Simple  calculation  of  this  angle  assuming  the  electrons 
move  at  relativistic  velocity  in  the  magnetic  field  established  by  the  vacuum 
beam  current,  indicates  an  entrance  angle  of  about  42°  to  the  normal.  This 
is  close  to  the  value  indicated  by  the  data  of  Figure  10,  assuming  major 
obscuration  over  about  3 cm  to  the  left  of  the  web  and  a web  depth  of  3 cm. 

With  dosimetry  film  placed  5 cm  away  from  the  foil  in  air,  these  large 
gradients  in  beam  current  measured  at  the  foil  were  not  detectable  due  to 
randomized  gas  scattering.  The  actual  experimental  conditions  called  for 
1.  7-3  atm  of  Argon,  which  further  reduced  the  gradients. 

Deposition  uniformity  in  the  beam  direction  into  the  high  pressure  gas 
was  perhaps  the  most  important  characteristic,  due  to  considerations  of 
discharge  stability  to  be  discussed  with  the  experimental  results.  In  order 
to  determine  the  extent  of  this  important  effect,  both  theoretical  calculations 
and  experimental  estimates  were  undertaken. 

The  theoretical  effort  was  based  on  an  electron,  Monte  Carlo  scattering 

computer  code  which  calculated  the  combined  effects  of  foil  scattering,  gas 

(34) 

scattering  and  secondary  emission  from  the  anode  . Results  for  a 350 
keV  beam  passing  through  a 50  pm  thick  foil  and  impinging  on  3 atm  of 
Argon  are  shown  in  Figure  11.  The  view  shown  is  made  to  simulate  the  end 
view  of  the  present  apparatus  which  had  a total  beam  aperture  of  20  cm  and  a 
maximum  anode-cathode  separation  of  10  cm.  The  contour  lines  indicate 


J.  A.  Halbleib,  Jr.  and  W.  H.  Vandevender,  "Cyltran:  A cylindrical- 
Gcometry  Multimaterial  Electron /Photon  Monte  Carlo  transport  Code", 
Sandia  Laboratories,  Albuquerque,  New  Mexico  87115. 
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points  of  constant  energy  deposition  in  keV/cm  which  when  multiplied 

2 3 

by  the  initial  current  density  in  A /cm,  indicate  power  deposition  in  kW/cm  . 

The  analysis  shows  that  a peak  in  the  deposition  occurs  at  about  4 cm  into  the 

gas  after  which  it  falls  off  by  25%  over  10  cm.  In  the  other  direction,  from 

center  to  edge,  the  fall  off  is  roughly  50%. 

An  approximate  experimental  measurement  of  beam  current  density 

fall-off  with  distance  from  the  aluminum  foil  was  attempted  by  placing  both 

dosimetry  film  and  a Faraday  cup  within  the  high  pressure  laser  cell.  The 

two  techniques  correlated  within  about  20%  and  indicated  a very  pronounced 

current  density  fall-off  as  shown  in  Figure  12.  The  beam  energy  was  320  keV, 

and  the  cell  was  pressurized  with  argon.  These  measurements,  of  course, 

indicate  current  density  only,  however,  an  independent  measurement  of  total 

energy  deposition  via  spatially  resolved  fluorescence  in  a separate  device  at 

(35) 

NRTC  yielded  substantially  similar  gradients  . 

It  will  be  shown  in  a subsequent  section  dealing  with  the  experimental 
results  how  nonuniformities  in  beam  energy  deposition  were  directly 
reflected  in  the  laser  output  energy  profiles  and  how  such  nonuniformities 
greatly  limited  the  attainable  discharge  energy  input. 

2.2  Laser  Plenum  and  Gas  Handling.  Design  of  the  laser  plenum  and  the 
associated  gas  handling  apparatus  was  guided  by  the  requirements  for  operation 
with  halogens  mixed  in  high  pressure  diluents.  Additionally,  the  plenum 
geometry  was  constrained  by  the  necessity  for  low  discharge  current  return 
path  inductance.  With  regard  to  materials,  it  was  considered  crucial  that 
those  appropriate  to  proper  passivation  and  opei'ational  cleanliness  be 
chosen.  Therefore,  all  metallic  parts  were  constructed  out  of  either 
stainless  steel  or  aluminum,  and  teflon  was  used  for  all  Insulating  areas.  In 

(35) 

J.  West,  Northrop  Research  and  Technology  Center,  private  communication. 
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order  to  minimize  discharge  circuit  inductance,  the  basic  plenum  design 
called  for  four  parallel  discharge  current  feedthroughs  and  a 300  liter  high 
pressure  cell  in  roughly  a semi-cylindrical  geometry. 

Operationally,  the  apparatus  performed  as  designed.  The  important 
limitation  to  the  geometry  was  foreseen  as  arcing  either  from  the  discharge 
electrode  to  the  chamber  walls  or  tracking  along  the  teflon  feedthrough 
insulators.  Neither  occurred  except  under  extreme  tests,  wholly 
unrepresentative  of  actual  experimental  conditions. 

Laser  Plenum.  The  high  pressure  laser  plenum  was  formed  from  a 
stainless  steel,  half-cylinder  attached  to  a massive  flange  which  provided  an 
interface  to  the  electron  gun.  Material  sizes  and  methods  of  construction 
were  consistent  with  a 10  atm  static  pressure  load,  assuming  fill  pressures 
up  to  5 atm  and  10  atm,  pulsed  over-pressures.  Safety  factors  were 
therefore  judged  to  be  a factor  of  five.  The  plenum  and  the  pumping  apparatus 
attached  to  it  were  in  turn  secured  to  an  oil  filled,  aluminum  container 
housing  the  discharge  switching  circuitry  and  capacitor  bank.  The  two 
attached  structures  were  then  allowed  to  roll  on  a track  providing  for 
separation  from  the  electron  gun  for  servicing  both  the  electron  gun  and  the 
plenum  interior.  These  components  are  shown  in  detail  in  Figures  13  and 
14  which  provide  an  end  and  front  view,  respectively.  A picture  of  the 
assembled  apparatus  is  shown  in  Figure  15. 

Figure  13  shows  the  method  of  attachment  of  the  plenum  to  the  discharge 

capacitor  bank  housing  as  well  as  the  placement  of  the  pumping  apparatus. 

With  regard  to  the  pumping  method,  a 6"  diameter  diffusion  pump  with  a liquid 

nitrogen  cold  trap  was  used  and  found  capable  of  pumping  the  entire  system 
-5 

down  to  5 x 10  mmHg  before  e '.ch  gas  fill. 

Figure  14  indicates  the  positions  of  the  four  discharge  current  feeds  and  a 
fifth  feed  which  allowed  direct  attachment  of  voltage  sampling  circuitry  to  the 
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Feedthrough 


Figure  15  - Picture  of  the  Device 


discharge  electrode.  Also  shown  is  the  diagonal  optical  path  for  the  He-Cd 
laser  used  to  monitor  partial  pressure  through  absorption  at  3250  A. 

These  diagnostics  are  discussed  below. 

Detailed  schematics  showing  the  relative  positions  of  the  important 
elements  are  shown  in  Figures  16  and  17.  Of  special  interest  are  the 
feedthrough  geometry,  anode  attachment  and  shape,  and  optical  attachment 
and  method  of  adjustment. 

Insulation  of  the  high  voltage  feedthroughs  inside  the  plenum  \v=’s  provided 

by  the  teflon  sleeving.  Outside  this  volume,  insulation  was  provided  by  the 

transformer  oil  which  was  allowed  to  penetrate  the  bellows  attachment  area 

from  the  discharge  capacitor  tank.  Attachment  to  the  discharge  anode  was 

made  via  a slip  collar,  allowing  for  adjustment  of  the  anode-cathode  dimension 

from  a maximum  of  10  cm  to  a minimum  of  2 cm.  The  anode  itself,  of 
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overall  dimension  20  x 100  cm  , was  formed  from  aluminum  into  a Rogowski 
profile  consistent  with  a 10  cm  electrode  spacing.  The  profile  geometry 
was  considered  successful  in  that  traces  of  arcing  at  all  anode-cathode 
separations  took  the  form  of  diffused  blemishes  along  the  central  third  of 
the  electrode,  and  no  tracking  was  observed  along  its  perimenter. 

Contributing  to  this  success  was  the  presence  of  a uniform  ground  plane 
presented  to  the  electrode  well  beyond  its  edges,  as  indicated. 

Placement  of  the  optical  extraction  area  in  relation  to  the  excitation 
volume  was  crucial,  owing  to  the  nonuniformity  of  energy  deposition  in  the 
high  pressure  gases.  As  indicated  in  the  accompanying  figures,  extraction 
was  obtained  immediately  beyond  the  ground  plane  screen,  however,  the 
extent  of  the  extraction  area  in  the  anode  direction  depended  upon  the  optical 
holder  profile.  In  this  regard,  two  types  of  holder  were  used  providing  for 
optical  path  lengths  of  1.2  and  2.0  M.  The  first,  shown  in  heavy  lines  in  the 


Figure  16.  Detailed  Plenum  Side  View 


Figure  17.  Detailed  Plenum  End  View 


figures,  provided  for  extension  of  a 10  x 20  cm  optic  well  inside  the  plenum 
volume  and  in  this  case,  the  available  distance  in  the  anode -cathode  direction 
open  for  extraction  was  8.4  cm.  The  second  method  used  a large,  28  cm 
diameter  optic  indicated  by  the  dotted  outline  in  Figure  16  attached  directly 
to  the  exit  port,  thus  providing  for  extrawcion  over  the  full  10  cm  anode-cathode 
Jis'  nee.  In  both  arrangements,  adjustment  of  the  optics  was  provided  by 
the  rnal  adjustment  bolts  attached  to  the  flexible  bellows  assembly.  The 
provisi..n  for  changing  optical  path  length  at  the  expense  of  total  extraction 
volume  was  important,  since  it  was  shown  experimentally  that  the  decrease 
in  path  length  indicated  above  enhanced  the  specific  energy  extraction  by 
about  50%. 

Gas  Handling.  The  gas  handling  system  was  designed  to  fulfill  several 
important  criteria.  The  entire  system  was  required  to  safely  handle  Fluorine 
at  high  pressures  in  the  storage  and  passivation  procedures.  This  also 
included  the  parallel  accomodation  of  high  pressure  rare  gases  without  con- 
tamination ard  with  additional  purification  if  indicated.  Proper  mixing  of 
the  various  gas  constituents  at  high  pressure  was  to  be  guaranteed  before 
each  electron  gun  bursf.  A continuous  monitor  of  the  F^^  partial  pressure 
was  to  be  provided  since  this  cj'ucial  constituent  was  present  in  small  amounts 
and  subject  to  chemical  reaction  at  the  chamber  surfaces.  These  various 
topics  are  described  below  in  detail  including  the  operational  experiences. 

It  was  found,  in  general,  that  the  entire  system  performed  well  above  the 
experimental  requirements. 

The  gas  storage  system  was  by  conventional  means  ’dth  an  explosion- 
proof  casing  for  the  Fluorine  and  standard  high  pressure  bottles  for  the  rare 
gases.  Commerically  available  gases  were  used  and  this  included  Argon 
(99.  998%),  Krypton  (99.995%),  Xenon  (99.  995%),  Neon  (99.  99%)  and 
Fluorine  (97%).  No  attempt  was  made  to  verify  the  analyses  provided  by 
supplier  and  Fluorine  was  not  purified  beyond  its  stated  level.  However, 


an  attempt  was  made  to  eliminate  water  vapor  and  hydrocarbon  contaminants 
through  the  use  of  alumina  gel  and  molecular  Tieve  filters.  It  was  found  that 
the  laser  behavior  did  not  change  with  the  addition  of  these  filter  types.  The 
gases  were  admitted  through  well  passivated  stainless  steel  tubes  to  the  laser 
plenum  where  mixing  took  place.  The  partial  pressures  of  Fluorine,  un  the 
order  of  2-4  Torr,  and  Krypton,  on  the  order  of  100  Torr,  were  measured 
with  a capacitance  minometer  which  used  a stainless  steel  diaphragm  in 
contact  with  the  gas.  The  higher  total  pressures,  on  the  order  of  2-4  atm, 
were  measured  with  stainless  s' oel  Bourdon  tube  gauges. 

A high  capacity  cryo-pump  was  used  to  reclaim  the  costly  rare  gases. 

The  pump  was  composed  simply  of  an  aluminum  cylinder  surrounded  r y a 
liquid  nitrogen  bath.  Gas  which  was  subjected  to  about  twenty  bursts  with  the 
electron  beam  and  reclaimed  in  this  manner  showed  no  degradation  due  to 
contamination.  Upon  re-use,  no  change  in  laser  performance  was  apparent. 

Using  a He-Cd  laser  source,  absorption  at  3250  k was  used  to  monitor 

partial  pressure  in  the  presence  of  high  pressure  rare  gases  in  the  laser 

plenum  before  each  electron  gun  burst.  The  monitor  also  provided  an 

indication  of  the  gas  mixing  time  after  high  velocity  injection  of  the  rare  gases 

through  the  small  inlet  orifice.  The  absorption  coefficient  of  0.  1 atm  ^ cm 

was  checked  for  the  1.  7 M absorption  path  length  with  the  capacitance  manometer 

(36) 

in  pure  Fluorine  at  low  pressures'  . For  these  low  concentrations,  total 
absorption  was  5% /Torr  of  F^,  well  within  the  resolution  of  the  apparatus. 

It  was  estimated  chat  accuracies  of  0.25  Torr  were  realised.  It  was  also 
determined  that  the  absorption  was  dependent  only  upon  the  F^  level  and 
independent  of  the  neutral  rare  gas  pressure.  Sines  the  technique  provided  an 
"on-line"  measurement  of  F^  partial  pressure,  disturbances  in  concentration 

J.  S.  Whittier,  M.  L.  Lundguist,  A.  Ching,  G.  E.  Thornton,  R.  Holland,  Jr., 
J.  Appl.  Phys.  £7,  3542  (1976),  and  private  communication. 
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along  the  absorption  path  were  discernible  during  the  nnixing  process.  Mixing 
was  considered  complete  when  the  absorption  level  measured  before  introduction 
of  the  rare  gases  was  attained.  Mixing  times  were  on  the  order  of  several 
minutes  and  operation  of  the  laser  at  times  greatly  in  excess  of  this  value 
showed  no  change  in  performance. 

The  arrangement  of  the  monitor  is  shown  in  Figure  18.  A 2 mW,  He-Cd 
laser  operating  at  3250  k was  used  as  the  cw  source  with  the  output  split  into 
two  paths,  one  propagating  in  air  as  the  reference  and  the  other  propagating 
diagonally  through  the  laser  plenum  adjacent  to  the  discharge  ground  screen. 

A chopping  wheel  was  used  to  allow  the  receiving  photodiode  to  sample  each 
beam  alternately.  The  photodiode  output  was  displayed  on  an  oscilloscope  and 
the  difference  signal  interpreted  directly  as  a measure  of  absorption. 

Absolute  calibration  of  the  system  was  easily  accomplished  with  the  plenum 
evacuated,  in  which  case  the  angle  sensitive  attenuator  was  adjusted  to  give 
zero  difference  signal. 


2.3  Discharge  Circuit.  Discharges  in  the  high  pressure  gas  mixtures 
typical  of  the  KrF  laser  are  characterized  by  very  low  resistance.  For 
devices  of  the  size  under  discussion,  the  characteristic  discharge  impedance, 
a function  of  time,  is  on  the  order  of  0.  25  ohm.  In  this  application, 
typical  high  voltage  pulsed  circuits  are  limited  by  their  inherent  inductance 
which  can  be  on  the  order  of  100-200  nH.  The  immediate  problem,  then,  is 
one  of  adapting  conventional  pulsed  discharge  circuitry,  which  normally 
requires  high  impedance  loads  for  efficient  energy  transfer.  Design  of  the 
present  discharge  cir  uit  was  therefore  guided  by  the  requirement  for 
minimization  of  total  circuit  inductance.  The  basic  circuit  itself  was  of  a 
simple  capacitive  discharge  type  with  subsequent  modification  into  a two  stage, 
pulse  charged  configuration  in  an  attempt  to  circumvent  a portion  of  the 
circuit  inductance. 
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An  extensive  circuit  analysis  was  first  undertaken  with  best  guess 
parameters.  Based  on  design  criteria  obtained  in  this  fashion,  the  actual 
circuit  was  constructed  and  tested  with  a fixed  dummy  load.  It  was  found 
that  the  circuit  behaved  substantially  as  predicted  and  that  the  total  circuit 
inductance  was  about  2Z5  nH  as  estimated. 

Circuit  Analysis.  The  discharge  circuit  model  is  shown  in  Figure  19 

with  a typical  set  of  component  values.  Three  parallel  sets  of  capacitive 

storage  are  indicated,  each  with  a series  inductance  of  50  nH  and 

resistance  of  0.  001  ohm.  The  load  is  represented  by  the  variable  resistance, 

R,  and  the  balance  of  the  circuit  including  the  switches  and  feedthroughs 

is  approximated  by  the  150  nH  inductance  and  0.  010  ohm  resistance.  The 

circuit  pictured  is  a Chebishev  filter,  modified  to  work  into  a low  impedance 
(37) 

load  . Basically,  the  various  capacitive  and  fixed  inductance  legs  can  be 
thought  of  as  contributing  different  frequency  components  to  the  final 
waveform  observed  at  the  load.  The  complication  arises  in  the  interaction 
of  the  various  legs  among  themselves.  Calculation  of  the  fairly  complex 
resultant  waveform  was  undertaken  with  a computer  code  which  correctly 
modeled  all  aspects  of  the  circuit.  The  component  values  were  extensively 
varied  theoretically  with  the  result  that  those  indicated  were  optimum. 

In  order  to  determine  the  circuit  sensitivity  to  inductance,  the  output  was 
calculated  for  various  experimentally  attainable  values  of  inductance; 
and  a load  resistance  of  0.5  ohm.  The  result  is  shown  in  Figure  20.  For 
an  initial  charge  voltage  of  100  kV,  a change  in  inductance  from  50  to  200 
nH  results  in  a factor  of  two  Increase  in  risetime  to  maximum  voltage,  and 
over  this  range,  the  maximum  achievable  voltage  applied  across  the 
discharge  drops  from  80  to  65%  of  the  initial  charge  voltage.  Considering 

(37) 

' Reference  Data  for  Radio  Engineers,  Fifth  Edition,  Howard  W.  Sams  & 
Co.,  Inc.,  p.  8-24  (1968). 
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a nominal  1 s pulse  length,  tne  total  input  energy  likewise  drops  from  8.  6 
to  5.  7 kJ.  Inductance  is  therefore  crucial  to  discharge  performance.  Figure 
21  shows  the  circuit  response  to  variable  load  impedance  assuming  a circuit 
inducatnce  of  150  nH  in  addition  to  the  50  nH  in  each  circuit  leg.  Since  the 
electron  beam  current  density  determines  the  load  impedance  through  primary 
ionization,  each  load  resistor  is  associated  with  a beam,  current  density. 
Additionally,  each  value  of  current  density  determines  the  pulse  length 
according  to  limitations  of  the  cold  cathode  electron  gun  technology.  For  the 
0.  5 ohm  curve,  the  2.  0 pF  capacitor  has  been  increased  to  4.  0 p F in  order 
to  minimize  voltage  droop  over  a 2 fi  s pulse.  The  combination  of  low  load 
impedance  with  a moderately  high  circuit  inductance  results  in  very  poor 
energy  coupling  into  the  discharge.  The  maximum  applied  voltage  is  seen 
to  drop  from  70  to  45%  of  the  charge  voltage  and  risetime  is  extended  from 
400  to  700  ns.  The  spans  shown  at  the  right  border  of  the  graph  indicate 
the  extent  of  optimum  E/N  for  2 and  3 atm  total  pressures  and  a 10  cm  anode- 
cathode  separation.  The  implication  is  that  even  in  the  presence  of  poor 
discharge  coupling,  the  experimental  parameters  are  within  a useful  range. 

An  attempt  was  made  to  circumvent  a portion  of  the  circuit  inductance  by 
placing  additional  capacitance  and  a series  switch  directly  across  the  load, 
as  indicated  in  Figure  22.  Here,  and  the  200  nH  series  inductance 
represent  the  initial  circuit,  and  its  50  nH  series  inductance  represent 
the  circuit  addition,  the  50  nH  inductor  in  series  with  the  load  represents 
the  remaining  circuit  inductance,  and  the  additional  switch  is  provided  by 
the  electron  beam  which  induces  an  instantaneous  impedance  collapse 
from  infinity  to  0.2  ohms  across  the  discharge.  Physically,  these  modifications 
were  carried  out  by  placing  capacitors  encased  in  Pyrex  cylinders  directly 
within  the  laser  plenum  beyond  the  major  current  feedthroughs  and  switches. 

The  discharge  was  operated  by  charging  the  main  bank,  C^,  switching  this 
onto  and  waiting  approximately  1 pi  s for  voltage  buildup,  and  then  firing 
the  electron  gun  which  allowed  current  to  flow  through  the  medium  from  the 
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composite  circuit.  The  circuit  was  analyzed  simply  by  assuming  both  C ^ 
and  charged  to  the  same  potential  and  switched  across  a fixed  load 
simultaneously.  Figure  23  shows  the  results  for  various  values  of  C^.  The 
indication  was  that  the  addition  of  greatly  improved  the  circuit  risetime 
and  that  the  value  of  should  be  as  large  as  possible.  Physical  constraints 
imposed  by  the  plenum  geometry  limited  to  0.  7 pF. 

Electrical  and  Mechanical  Arrangement.  The  discharge  circuit  design 
was  guided  primarily  by  the  need  for  low  inductance  and  secondarily  by  the 
requirement  for  reliability  and  versatility.  Therefore,  the  basic  layout 
called  for  stripline  capacitive  storage,  low  inductance  switching  techniques 
and  multiple  current  feedthroughs.  Additionally,  the  design  required  that  the 
circuit  components  be  kept  in  close  proximity,  therefore,  the  entire  circuit 
was  submerged  in  insulating  oil.  The  relative  arrangement  of  the  circuit 
components  is  shown  in  Figure  24. 

The  capacitor  ba.kk  was  composed  of  twenty  parallel  capacitors  in  five 
rows  connected  by  metal  sheets  in  a stripline  configuration.  The  row  closest 
to  the  main  switch  contained  four  capacitors,  two  each  of  0.015  and  0.  15  pF 
capacitance  and  the  remaining  pieces  were  each  0.25  pF.  The  intent  was  to 
minimize  inductance  between  the  smallest  capacitances  and  the  load  and 
possibly  add  an  inductor  to  the  remaining  bank,  to  slow  up  that  portion  of  the 
circuit  and  thus  keep  the  voltage  from  drooping,  in  the  spirit  of  the  Chebyshev 
circuit  described  above.  However,  under  operation  it  was  found  that  significant 
voltage  droop  was  desirable  in  order  to  keep  the  discharge  from  arcing, 
therefore  this  modification  was  not  attempted. 

The  switching  circuitry  was  required  to  meet  several  criteria:  jitter 
after  application  of  the  trigger  pulse  was  to  be  under  20  ns,  reliable  switching 
was  to  be  provided  for  initial  voltages  from  20-100  kV  and  total  currents  of 
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Figure  23.  Pulse  Charged  Circuit  Response 
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Figure  24.  Discharge  Circuit  Arrangement 


100  kA,  and  switch  inductance  was  limited  to  20  nH,  In  order  to  meet  these 
specifications,  the  most  reasonable  approach  was  a rail  gap  composed  of  two 
continuous  electrodes.  Identical  geometries  were  used  for  both  the  main  and 
crowbar  switches.  The  other  alternative  of  individual  spark  gaps  in  parallel  with 
transit  time  isolation  was  considered  far  less  versatile  and  more  difficult 
to  accomodate  physically. 

Commercially  available  trigger  generators  and  rail  gaps,  supplied  by  the 
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Maxwell  Laboratories  were  used'  . The  rail  gaps  were  constructed  in  a 30  cm 
long  three  electrode  geometry  with  a sharp-edged  continuous  trigger  electrode 
placed  at  the  1/3-2/3  point  between  two  cylindrical  main  electrodes.  Two 
electrode  spacings  were  required  corresponding  to  the  voltage  ranges 
20-50  kV  and  50-100  kV,  and  the  gaps  operated  in  a 1/5.7  mixture  of  SF^ 
in  Argon  at  pressures  from  20-60  psig.  The  trigger  generators  produced 
a 60  kV  peak  voltage  pulse  at  a differential  risetime  of  10  kV/ns  and  were 
sufficient  to  cause  the  formation  of  multiple  current  channels  ( > 6)  over  the 
full  electrode  extent  for  the  entire  voltage  range.  The  trigger  pulses  were 
capacitively  coupled  into  each  of  the  resistively  biased  trigger  electrodes. 

In  addition,  an  isolating  jpark  gap  and  capacitive  bias  in  the  crowbar  trigger 
input  guaranteed  isolation  and  reliable  switch  initiation  after  significant 
voltage  droop  late  in  the  pulse.  It  was  found  that  with  this  arrangement,  the 
crowbar  switch  would  reliably  hold  off  100  kV  at  the  beginning  of  a pulse 
and  reliably  '.rigger  after  the  voltage  had  drooped  to  30  kV,  providing  for  a 
wide  margin  of  operation  at  other  voltages  as  well.  Jitter  for  both  the  main 
and  crowbar  switches  was  measured  at  under  20  ns  over  the  parameter  ra.iges 
indicated  above. 


As  outlined  in  the  section  dealing  with  plenum  design,  energy  from  the 
capacitive  storage  bank  was  transmitted  to  the  discharge  electrode  via  fou: 


(38)  , 


Maxwell  Laboratories,  San  Diego,  California. 
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high  voltage  feedthroughs.  Briefly,  the  feedthrough  length  was  60  cm  and 
coaxial  over  about  20  cm  with  the  remainder  open  to  either  the  high  pressure 
laser  gases  or  insulating  oil. 

Various  peripheral  circuit  elements  were  added  for  convenience.  A 10  ohm 
Ammonium  Chloride  liquid  resistor  was  connected  at  the  output  of  the  main 
switch  to  ground  in  order  to  provide  periodic  calibration  of  the  voltage  and 
current  monit'  is  and  'milage  bleed-off  in  the  event  of  false  triggering 
without  electron  beam  operation.  A mechanical  crowbar  was  also  attached 
at  the  same  circuit  point  to  prevent  premature  voltage  application  to  the 
discharge  electrode. 


Basic  circuit  parameters  were  measured  with  a voltage  probe  attached 

directly  to  the  discharge  electrode  and  current  monitor  in  the  ground  leg  of  the 

capacitor  bank.  Both  monitors  were  designed  to  provide  high  level  signals 

without  direct  contact,  thereby  eliminating  sources  for  interference.  The 

voltage  probe  consisted  of  a 500  ohm  fixed  resistor  immersed  in  insulating 

oil  within  a teflon  cylinder.  The  cylinder  was  fit  inside  the  plenum  and  the 

resistor  attached  through  a conducting  cylinder  end  cap  to  the  discharge 

electrode.  The  opposite  end  of  the  resistor,  likewise  connected  to  an 

eicternal  cylinder  end  cap,  was  connected  to  ground  thiough  a toroidal 
(33) 

Pearson  ' current  probe  with  a risetime  of  20  ns.  Since  the  Pearson  probe 
calibration  was  IV/A,  signal  levels  were  on  the  order  of  50-100  V.  Absolute 
calibration  and  frequency  response  were  verified  by  parallel  measurement 
with  a Tektronix  voltage  probe  and  the  system  operating  at  reduced  voltage. 


The  current  monitor,  provided  by  Physic:.  International^"*^^,  was  basically 
a Rogowski  coil  wound  on  a line  1/2  M long  and  placed  within  a notch  in  the 
condc  ting  sheet  from  the  capacitor  termuial  to  ground,  thus  providing  for 
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Physics  International,  San  Diego.  -Talifcrnia. 
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a low  inductance  connection.  The  linear  current  probe  responded  to  the 
derivative  of  magnetic  field  and  required  an  integrator  to  display  current  versus 
time.  The  system  had  a risetime  of  10  ns  with  an  integrator  of  time  constant 
20  fJi  s (time  to  droop  10%).  Calibration  of  the  system  was  checked  with  a 
reduced  capacitor  bank  and  the  output  measured  with  a Pearson  probe  through 
a fixed  dummy  load  resistor. 

Performance.  The  entire  discharge  supply  was  tested  with  a 1 ohm 
NH^Cl  liquid  resistor  as  a dummy  load.  In  order  to  determine  the  inductance 
under  actual  operating  conditions,  the  load  was  formed  from  four  parallel  legs 
7.  5 cm  in  diameter  by  30  cm  long,  and  returned  to  ground  with  IM  wide 
sheets  to  simulate  the  plenum  walls.  Typical  voltage  and  ;urrent  waveforms 
are  shown  in  Figure  25.  In  Figure  25a,  the  trace  starting  at  the  upper  left  and 
going  negative  is  total  current  measured  by  the  linear  probe,  and  the  bottom, 
positive  going  trace  is  voltage  applied  to  the  dummy  load  as  measured  by 
the  voltage  monitor.  In  this  case,  bank  capacitance  was  2.  33  pF  with  an 
initial  charge  voltage  of  50  kV,  and  the  crowbar  was  set  to  actuate  3 p s after 
pulse  initiation.  The  crowbar  action  is  clearly  apparent  at  the  point  where 
the  total  current  begins  to  rise  steeply  and  load  voltage  begins  to  fall  abruptly. 
With  the  traces  extended  in  time  as  in  Figure  25b,  an  indication  of  the 
ringing  frequency  through  the  crowbar  circuit  was  obtained.  The  crowbar 
inductance,  calculated  from  the  usual  expression,  was  about  174  nH.  In 
similar  fashion,  the  total  discharge  circuit  inductance  was  measured  to  be  about 
225  nH.  Therefore,  it  was  expected  that  the  actual  discharge  with  impedance 
less  than  0.  2 Q would  compete  favori.bly  with  the  crowbar  for  capacitor  bank 
energy  and  indeed  this  was  the  case.  Discharge  arcing  was  found  to  take 
place  via  a number  cf  low  inductance  diffuse  glows  along  the  discharge 
electrode,  and  did  not  enhance  possible  diversion  of  energy  through  the 
crowbar.  Therefore,  although  applicable  to  some  situations,  the  crowbar 
circuit  was  generally  not  useful. 
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Figure  25.  Capacitive  Discharge  Circuit  Waveforms 


Typical  traces  for  the  pulse  charged  circuit  operating  into  a plasma 
generated  by  the  electron  beam  are  shown  in  Figure  26.  In  this  case, 

C = 2.  33  |JL  F,  = 0.  7 |i  F and  the  initial  charge  voltage  was  30  kV. 

Figure  26a  shows  the  electron  gun  voltage  which  approximates  beam  current 
also,  and  Figure  b shows  the  corresponding  voltage  (upper  trace)  and  current 
(lower  trace).  The  discharge  electrode  voltage  is  characterized  by  a 
steeply  rising  leading  edge  as  capacitor  C^,  initially  at  ground  potential, 
begins  to  charge.  Likewise,  the  current  waveform  indicates  this  charging 
cycle  over  the  first  600  ns.  Abruptly  at  600  ns  into  the  pulse  charge,  the 
electron  beam  current  is  sufficient  to  cause  a drastic  drop  in  the  discharge 
impedance  from  inf  nity  to  about  0.5  ohm  at  which  point  the  discharge  voltage 
drops  to  its  operating  point,  at  10  kV.  After  this  point,  the  current  probe, 
which  measures  total  current  in  the  ground  leg  of  the  main  capacitive  storage, 
does  not  sense  true  discharge  current  supplied  by  both  and  Discharge 
voltage,  on  the  other  hand,  rises  then  falls  as  the  main  bank  is  depleted. 
Finally,  at  termination  of  the  electron  gun  pulse,  the  discharge  impedance 
rises  and  this  is  indicated  in  the  voltage  trace  as  a slight  peak  at  1.4  p s. 

2.  4 Laser  Diagnostics  and  Optics.  The  laser  output  beam  was  diagnosed 
in  its  important  aspects  through  a variety  of  techniques.  Total  beam  energy 
was  measured  calorimetrically,  and  two  types  of  calorimeter  were  tested  in 
order  to  provide  a consistency  check.  The  beam  temporal  behavior  was 
monitored  by  a fast  response  photodiode  with  ultraviolet  filters  at  the  input  to 
shield  against  extraneous  flourescence.  Finally,  approximate  but  effective 
beam  profile  measurements  were  obtained  through  beam  burn  patterns  on 
unexposed  and  developed  polaroid  film.  These  techniques  proved  highly 
consistent  and  reliable. 

The  laser  optics  involved  state  of  the  art  technology  to  produce  high 
damage  threshold,  low  absorptance  reflectors.  Since  periodic  optical  damage 
was  caused  by  the  combined  effects  of  intense  laser  radiation,  F^  attack  and 
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Figure  26.  Pulse  Charged  Circuit  Waveforms 
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x-ray  radiation,  optical  quality  was  continually  monitored  by  a simple, 

2500  A reflectometer.  In  order  to  reduce  costs  incurred  as  a result  of 
large  area  optic  damage,  a set  of  small  area  optics  was  used  in  the  initial 
parameterization  studies.  Asa  result  of  these  precautions  and  careful 
evaluation  of  various  vendors,  performance  of  the  optics  was  entirely 
satisfactory. 

Electronics.  Design  of  the  electronics  arrangement  was  guided  primarily 
by  the  need  to  eliminate  emi  from  both  the  triggering  circuits  and  diagnostic 
electronics.  Primary  shielding  of  the  diagnostic  electronics  was  provided 
by  a solid  wall  shielded  room  with  100  db  attenuation  at  100  MHz.  All  of  the 
oscilloscopes,  photodiode  bias  supplies  and  low  level  delay  trigger 
circuitry  were  placed  in  this  enclosure.  The  various  diagnostic  elements 
themselves  were  located  at  the  appropriate  positions  external  to  the  shielded 
room  and  connected  to  it  through  grounded,  solid  conduit  and  shielded 
cable.  The  generation  of  false  trigger  signals  due  to  emi  in  the  various 
trigger  generators  was  a major  concern.  Therefore,  an  effort  was  made  to 
maintain  high  level  signals  at  the  components  in  the  emi  environment  and  to 
maintain  electrical  isolation  of  the  trigger  systems.  These  efforts 
proved  highly  successful,  since  false  triggering  was  not  experienced  and 
noise  levels  in  the  shielded  room  were  kept  to  approximately  two 
oscilloscope  trace  widths. 

Signals  were  recorded  on  100  MHs  bandwidth  Tektronix  oscilloscopes, 
and  common  external  triggering  was  provided  by  a synchronizing  pulse 
from  the  electron  gun  trigger  generator.  Coincidence  of  the  various 
oscilloscope  traces  was  within  20  ns  and  this  was  considered  adequate 
at  the  200  ns/div  writing  speeds  used. 
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Diagnostics.  All  the  laser  diagnostics  were  arranged  on  a massive 
optical  table  as  shown  in  Figure  27.  The  laser  beam  at  full  power  was 
transported  first  to  a large  uncoated  beam  splitter  which  diverted  10%  of  the 
power  into  the  diagnostic  setup  and  allowed  the  balance  to  impinge  on  the  burn 
pattern  polaroid  film.  The  low  level  beam  was  then  partially  focused  by  an 
uncoated,  IM  focal  length  lens  onto  the  calorimeter,  either  a Scientech  #364 
or  Gentec  #ED  500.  Again,  10%  of  the  beam  entering  the  calorimeter  was 
diverted  into  a photodiode,  ITT  #4000(55).  The  10  cm  diameter  aperture  of 
the  Scientech  calorimeter  provided  for  complete  capture  of  the  partially 
focused  beam;  and  surface  absorptance  at  2500  A was  94%.  The  diode  was 
operated  at  a bias  of  2.5  kV  and  sufficient  attenuation  was  used  to  keep  the 
current  well  below  saturation.  In  the  absence  of  lasing  or  with  the  laser 
beam  blocked  from  entering  the  photodiode,  no  signal  was  detected, 
indicating  that  neither  fluorescence  nor  extraneous  signals  contributed  to 
the  observations.  All  of  the  beam  diagnostic  optics  were  made  of  high  purity 
quartz  with  absorption  less  than  0.  1%  in  the  ultraviolet  region  for  wavelengths 
greater  than  2300  k. 

The  beam  splitting  and  focusing  arrangement  was  calibrated  by  operating 
the  laser  at  greatly  reduced  power  and  allowing  the  full  energy  to  enter  the 
calorimeters  directly.  Additionally,  the  calorimeters  were  found  to  give 
identical  results  per  their  stated  calibrations  under  the  same  output 
conditions  as  indicated  by  the  photodiode. 

Optics.  The  laser  optical  system  was  formed  from  three  geometrically 
distinct  sets  of  optics.  (1)  Small  diameter  optics  for  the  initial  laser 
parameterization  studies.  These  were  7.  6 cm  in  diameter  allowing  for  a 
clear  aperture  of  25  cm^;  (2)  Large  area  rectangular  optics  for  full  scale 
extraction  with  the  shortest  possible  optical  cavity  length  of  1.2  M.  These 
provided  for  a clear  aperture  of  9 x 14  cm^  although  partially  obscured  by 
the  gun  foil  support  as  described  in  that  section;  (3)  Large  area  circular 


56 


optics  for  full  scale  extraction  at  the  largest  aperture  possible  but  with  an 
excessive  optical  cavity  length  of  1.5  M.  Total  optic  diameter  was  28  cm 
and  the  extraction  area  was  sufficient  to  sample  the  full  20  liter  gain  volume. 

In  all  cases,  the  optical  components  were  formed  from  high  purity  quartz  with 
multi-layer  dielectric  coatings  of  MgF^  and  ThF.  In  some  cases,  anti- 
reflection  coated  windows  of  ♦he  same  materials  were  used,  but  the  usual 
arrangement  was  with  reflectors  only  at  normal  incidence.  The  mirror 
surfaces  were  either  flat  or  possessed  a 10  M radius  of  curvature  to  provide 
for  slightly  stable  optical  resonators.  It  was  found  that  output  energy  and 
temporal  behavior  did  not  change  with  either  flat/flat  or  flat/radius 
resonators  and  since  beam  divergence  was  not  an  important  issue,  both  types 
of  resonator  were  used  interchangeably.  The  outward  force  of  the  high  pressure 
gas  induced  added  curvature  to  the  larger  area  optics  producing  radii  of 
curvature  estimated  at  20  M,  and  this  effect  may  have  contributed  to  the 
insensitivity  of  the  output  on  resonator  design. 

The  mirror  surfaces  were  susceptible  to  attack  by  F^  and  changed  character 
through  extended  use.  As  a result  it  was  considered  imperative  that 
reflectivities  be  checked  periodically,  approximately  every  100  shots.  Since 
the  large  size  of  the  mirrors  precluded  the  use  of  commercially  available 
instruments,  a simple  spot-check  (2500  + 50  k)  reflectometer  was  set  up  as 
shown  in  Figure  28.  The  intensities  of  the  sample  and  reference  beams  were 
obtained  by  alternately  blocking  them.  The  reflectivity  of  the  test  piece  was 
then  given  by  (1-R)r/R,  where  R is  reflectivity  of  the  partial  reflector  and  r 
is  the  ratio  of  sample  beam  to  reference  beam  intensities.  The  use  of  an 
aperture  over  the  ultraviolet  source  allowed  spatial  scan  of  the  optic  under 
test.  Calibration  tests  with  small  diameter  optics,  which  were  checked  in 
commercial  instruments,  indicated  accuracies  of  + 5%  in  reflectivity. 

Alignment  of  the  laser  resonator  under  operating  conditions  was  performed 
with  a spatially  filtered  and  collimated  He-Ne  laser  beam.  The  reflections  from 
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i.*ie  laser  mirror  interior  surfaces  were  propagated  over  approximately 
60  M with  resultant  spot  diameters  of  about  0.  5 cm.  Alignment  accuracy  was 
estimated  at  about  half  a spot  diameter,  or  8 sec  of  arc.  This  accuracy  was 
considered  to  be  well  within  acceptable  limits  since  degradation  of  the  laser 
output  required  resonator  misalignment  of  at  least  30  sec  of  arc. 

3.  J^ASER  KINETIC  ANALYSIS  AND  MODELING 

In  order  to  optimize  the  performance  of  theKrF  laser  and  to  understand 
its  scaling  to  high  power  and  high  volume,  it  is  necessary  to  understand  the 
kinetic  and  pumping  processes  which  are  responsible  for  the  formation  and 
quenching  of  the  excited  species.  In  addition,  it  is  necessary  to  analyze 
discharge  stability  and  effects  of  nonuniformity.  The  details  of  an  extensive 
analytical  model  and  comparisons  of  its  predictions  with  experimental  data 
from  the  large  scale-up  device  discussed  above  will  be  described. 

The  theoretical  model  for  an  electrically  excited  laser  requires  a 
coupled  analysis  wnich  includes  electron  kinetics,  molecular  kinetics, 
external  driving  circuit,  and  optical  radiative  extraction.  The  components  of 
such  an  analysis,  labelled  with  brief  descriptive  remarks,  are  summarized 
in  Figure  29.  A comprehensive  computer  code,  which  has  been  developed 
to  implement  this  coupled  analysis  for  a broad  class  of  electrically  excited 
laser  systems,  will  be  described  below.  The  molecular  kinetics  are 
described  by  the  master  equation  for  the  population  densities  of  all  the  atomic 
and  molecular  species  present  in  the  electrically  excited  gas  mixture,  in- 
cluding electrons,  ions,  neutrals,  excited  states,  molecules,  excimers,  etc. 

The  electron  kinetics  are  described  by  the  Boltzmann  transport  equation  for 
the  electron  energy  distribution  function  f(u),  from  which  all  plasma  parameters, 
secondary  electron  excitation  rates,  power  partitioning,  and  discharge 
efficiency  can  be  obtained.  The  formulation  of  the  Boltzmann  equation 
includes  elastic  and  inelastic  electron-neutral  collisions,  momentum 

'^This  analysis  is  formulated  for  a spatially  homogeneous  medium,  with  a 
simplified  model  of  the  optical  resonator,  which  is  assumed  to  be  formed  by 
two  plane  parallel  mirrors. 
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• MOLECULAR  KINETICS 

Master  Equation  for  population  densities  of  all  atomic  and  molecular 
species  (electrons,  ions,  neutrals,  excited  states,  excimers,  etc. ). 

• ELECTRON  KINETICS 

Boltzmann  transport  equation  for  the  electron  energy  distribution 
to  determine  plasma  parameters,  electrical  excitation  rates,  power 
partitioning,  and  discharge  efficiency. 

• CIRCUIT  ANALYSIS 

External  circuit  equations  to  describe  the  electrical  power  loading  In 
a low  impedance  gas  discharge,  whose  conductivity  is  function  of  t. 

• OPTICAL  EX'RACTiON 

Oscillator  analysis,  formulated  in  terms  of  intracavity  photon  density, 
with  transient  buiid-up  of  laser  mode  from  spontaneous  emission  and 
subsequent  relaxation  to  quasi-steady  condition  with  gain  at  threshold. 


Figure  29.  Coupled  Laser  Kinetics  Analysis 
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transfer  with  recoil,  electron-electron  scattering,  and  properly  incorporates 
superelastic  collisions  from  excited  species.  External  circuit  equations  are 
required  to  correctly  describe  the  electrical  power  loading  from  an  RLiC 
driving  circuit  into  a low  impedance  gas  discharge  whose  conductivity  is  a 
function  of  time.  Finally,  radiative  extraction  is  formulated  in  terms  of  the 
photon  density  in  the  cavity,  and  includes  a source  term  from  spontaneous 
emission,  amplification  from  stimulated  emission,  and  losses  from  output 
coupling  and/or  intracavity  absorption.  Details  of  each  component  of  the 
coupled  analysis  will  be  discussed  in  more  detail  below. 


3.  1 Molecular  Kinetics.  For  a given  set  cf  reactions,  the  master 
equation  for  the  population  densities  of  all  interacting  species  can  be 
constructed  as  shown  in  Figure  30.  For  any  reaction  involving  species 
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where  k,  and  k are  forward  and  reverse  rate  constants,  and  H.,  v.  are 
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(nonnegative  integer)  coefficients,  the  net  rate  is  given  by 
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where  [X*]  is  the  population  density  (cm  ) of  the  i species  X^.  The 
contribution  that  such  a reaction  makes  to  the  production  (or  loss)  of  the  i 
species  is  given  by 


d/dt  [X.]  = (v.  - p.)  R (3) 

1 1 1 

Thus,  for  an  arbitrary  reaction  scheme  containing  several  collision  processes 
(labeled  by  or),  the  complete  master  equation  becomes 
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(4) 


d/dt  [X.]  = ^ [v.(of)  - (of)  ] R (a) 

O'  ^ ^ 

It  is  impractical  to  construct  new  computer  codes  for  the  analysis  of  each  new 
kinetic  reaction  scheme.  The  straightforward  approach  to  write  the  subroutines 
which  define  the  molecular  kinetics  for  a complicated  reaction  scheme  is,  in 
general,  a difficult  and  time-consuming  task,  and  the  resulting  computer 
program  would  have  little  flexibility  for  analysis  of  any  system  except  those 
in  a very  limited  class.  (For  example,  there  are  currently  about  80  reactions 
considered  important  for  analysis  of  the  KrF  laser.  ) Furthermore,  as 
understanding  of  the  reaction  scheme  evolves,  addition  of  new  reactions  (or 
deletion  of  old  ones)  may  be  continually  required,  in  addition  to  simple 
modification  of  rate  constants  for  the  reactions  retained.  Such  a code  would 
itself  have  to  be  revised  in  a continuing  manner,  rather  than  merely  executed 
with  revised  values  for  the  rate  c instants. 

Therefore,  in  order  to  have  a powerful  analytical  capability  which  is 
applicable  to  a broad  class  of  problems,  a generalized  computer  cods  which 
itself  synthesizes  a coupled  analysis  (as  described  in  Figure  29)  has  been 
developed.  This  code  c utomatically  generates  its  own  subroutines  for  analysis 
of  the  molecular  kinetics  by  translating  symbolic  reactions  into  compu.sr- 
coded  equations.  Thus,  for  the  most  complicated  reaction  scheme  containing 
an  arbitrary  number  of  kinetic  collision  processes  and  interacting  species,  it 
is  possible  to  construct  the  complete  computer  code  required  with  virtually 
no  effort.  In  order  to  generate  the  ^oupled  analysis  for  any  reaction  scheme, 
all  that  is  required  is  to  provide  an  input  deck,  arbitrarily  long, 
consisting  of  pairs  of  cards,  the  first  containing  a reaction,  followed  by  a 
second  containing  forward  and  reverse  rate  constants  and  a reference.  The 
syntax  for  the  reactions  js  very  flexible,  with  a free  format  that  specifies 
the  reaction  just  as  it  would  normally  be  written.  The  content  of  each  reaction 
is  analyzed,  and  the  appearance  of  each  new  species  is  recognizea  and  its 
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name  stored.  The  syntax  of  each  reaction  is  subjected  to  numerous  tests  to 
detect  errors.  If  the  reaction  is  determined  to  be  acceptable,  it  is  translated 
into  appropriate  expressions  in  the  generation  of  synthesized  subroutines. 
Otherwise,  diagnostic  comments  are  produced.  Program  generation  and 
execution  are  protected  by  automatic  exit  (by  input  request)  if  specified  error 
conditions  are  encountered. 

Figure  31  presents  a schematic  flow  diagram  of  the  present  approach. 

The  rate  constants  initially  assumed  in  the  generation  of  the  program  can  be 
changed  in  the  subsequent  execution,  if  desired.  The  main  purpose  of  the 
rate  constants  in  the  initial  input  deck  is  to  define  whether  the  forward 
and/or  reverse  process  is  to  be  included  in  the  rate  expressions,  for  if  a 
zero  rate  constant  is  encountered  for  any  process  (other  than  secondary 
electron  collisions),  no  translation  of  the  forward  (or  backward)  term  occurs 
in  the  generated  subroutines.  If  the  correct  values  of  rate  constants  are 
known,  they  may  be  entered  for  once  and  for  all  in  the  original  reaction 
input  deck.  As  Figure  31  indicates,  the  synthesized  molecular  kinetics 
subroutines  are  automatically  combined  with  the  master  executive  program 
and  all  other  required  subroutines  to  form  a completely  self-contained 
coupled  analysis  based  upon  the  specified  reaction  scheme.  The  section 
enclosed  in  the  upper  box  in  Figure  31  is  thus  required  only  when  a new  reaction 
scheme  is  introduced.  For  all  subsequent  calculations  based  on  the  given 
reaction  scheme,  the  initial  conditions,  experimental  parameters,  rate 
modifications  (if  any),  and  I/O  and  control  parameters  are  entered,  and  the 
complete  code  constructed  in  the  upper  box  is  executed.  That  is,  only  the 
lower  section  of  Figure  31  is  involved  thereafter. 

There  are  several  obvious  advantages  to  this  approach.  First  of  all,  there 
is  the  simplicity  of  automation  and  the  minimization  of  the  possibility  of  error. 
For  example,  approximately  80  reactions  are  currently  believed  to  contribute 
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Synthesis  ; 


I 


Figure  31.  Generalized  Synthesis  Laser  Kinetics  Code 


to  the  KrF  reaction  kinetics,  and  to  write  the  subroutines  for  such  a 
complicated  system  would  not  be  simple.  Secondly,  the  program  diagnoses 
error  conditions  in  the  reaction  syntax  which  may  not  have  been  noticeable 
or  which  may  have  been  overlooked.  For  example,  duplicate  reactions  are 
detected  (even  when  written  backwards),  charge  particle  and  heavy  particle 
conservation  is  insured,  detail  balance  relations  for  binary  collision  processes 
are  enforced,  and  miscellaneous  other  error  conditions  are  detected. 

Secondary  electron  collision  processes  are  recognized  and  properly  coupled 
to  the  electron  kinetics  analysis,  except  for  the  case  with  no  discharge,  in 
which  case  their  rate  constants  can  be  defined  by  input.  Finally,  the 
subroutines  are  constructed  in  such  a way  that  null  operations  are  completely 
eliminated  (if  no  rate  constant  is  provided  for  the  forward  or  reverse  process, 
no  translation  of  the  corresponding  term  occurs),  and  in  such  a way  as  to 
optimize  execution  efficiency  (repetitious  or  unnecessary  multiplications 
are  minimized). 

It  is  usually  the  case  that  the  rate  constants  involved  in  molecular  kinetic 
reaction  schemes  can  vary  over  several  orders  of  magnitude,  and  therefore, 
the  resulting  master  equations  often  become  a "stiff"  system  of  differential 
equations.  Therefore,  the  approach  which  has  been  taken  is  to  employ  a 
multistep  integration  technique  developed  by  Gear^^*^^  for  solution  of  such 
equations.  This  method  automatically  adjusts  the  integration  step  size  as 
the  solution  proceeds,  in  such  a way  that  required  accuracy  conditions  are 
maintained.  The  Gear  method  requires  subroutines  not  only  for  the  rates  of 
change  of  the  population  densities,  dn./dt,  but  also  for  the  Jacobian, 

5n  /6n  , as  a function  of  time.  Because  numerical  evaluation  of  the  Jacobian 

1 j 

is  not  always  satisfactory,  it  is  necessary  to  generate  both  such  subroutines 
symbolically  in  the  synthesis  section  of  the  program,  where  the  reaction 

(40) 

C.  William  Gear,  Numerical  Initial  Value  Problems  in  Ordinary  Differential 

Equations.  Prentice-Hall,  Englewood  Cliffs,  N.  J.  (1971).  (Cf.  pp.  155-167). 
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scheme  is  translated  into  computer-coded  equations. 


The  entry  of  control  parameters,  experimental  parameters,  revised  rate 
constants,  initial  conditions,  output  requests,  etc.,  is  quite  flexible, 
and  permits  the  code  to  be  executed  for  a variety  of  situations  of  interest. 

In  the  case  of  pure  e-beam  excitation,  rate  constants  for  secondary  electron 
processes  (which  are  normally  coupled  to  the  Boltzmann  anai/sis)  default 
automatically  to  zero,  but  can  be  specified  by  input,  if  desired.  (This  is 
useful  for  electron  recombination  or  attachment  processes,  for  example.  ) 

This  makes  it  possible  to  use  the  same  general  code  for  both  discharge  and 
e-beam  excitation  conditions.  The  integration  of  the  coupled  set  of  equations 
over  the  total  specified  pulse  length  is  carried  out  with  the  Gear  technique 
which  automatically  adjusts  its  step  size.  However,  the  total  pulse  length 
is  divided  into  50  subintervals,  at  which  times  the  electron  kinetics  are 
updated,  and  a variety  of  output  options  can  be  specified. 

3.2  Radiative  Extraction.  Basic  features  of  the  radiation  analysis  are 

summarized  in  Figure  32.  It  is  convenient  to  be  able  to  treat  stimulated 

emission  processes  as  well  as  absorption  of  the  laser  radiation  by  species  in 

the  active  medium,  in  a manner  parallel  to  that  for  kinetic  collision  reactions. 

Therefore,  the  equation  for  radiative  extraction  is  formulated  in  terms  of 

the  photon  density  in  the  cavity  (n  , = I/chv),  and  includes  a source  term 

Ph 

from  spontaneous  emission,  amplification  from  stimulated  emission,  and 
losses  from  output  coupling  and/or  intracavity  absorption: 

n = (n^  ) (^/4tt)  + (L  /L  )c  [or  (t)  - a ] n . (5 

ph  IsrF  sp  g c th-"  ph 

and  are  the  lengths  of  the  gain  medium  and  optical  cavity,  = 

[y  + 1/2  ln(l/R)]/L  is  the  threshold  loss  coefficient  expressed  in  terms  of 

s 

the  intracavity  loss  per  pass  y and  the  output  coupling  reflectivity  R,  (fl/4n) 
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• Spatially  uniform  intracavity  optical  fields. 

• Plane  parallel  optical  resonator  (100%,  R),  separation  L^. 

• Spatially  uniform  gain  medium  of  length  Lg. 

• Distributed  output  coupling  and  loss  (Y)  coefficients. 

• Cavity  threshhold  gain  coefficient-. 

• [Y  (1/2)  ln{l/R)]/Lg 

• Buildup  of  optical  fields: 

Spontaneous  emission  in  solid  angle  0 defined  by  output 
coupling  aperture  and  resonator  length.- 

Amplification  from  stimulated  emission,  offset  by  the  in- 
tracavity absorption  and  output  coupling  losses. 

• Formulation  in  terms  of  intracavity  photon  density.- 


dHp^ltl/dt  • (Lg/L,lcla(t) 


Cavity  buildup  time:  t^.  • L^/Lga([.,c 


Figure  32.  Optical  Extraction 


is  the  fractional  part  of  the  total  spontaneous  emission  which  builds  up  the 
laser  mode  (taken  to  be  the  solid  angle  *1  subtended  by  the  output  coupling 
mirror  viewed  from  the  opposite  end  of  the  cavity),  and  or(t)  is  the  instantaneous 
net  gain  coefficient  in  the  medium  (i.  e. , gain  offset  by  absorption). 


o,(t)  = - I „^(t) 

A 


(6) 


where  the  sum  is  over  all  absorbing  species  A in  the  medium.  This  approach 
allows  the  radiation  intensity  to  be  described  for  situations  where  pumping 
and  kinetic  times  are  comparable  to  the  cavity  decay  time,  t^  = (L^ /L^or^j^c). 
For  such  situations,  gain  relaxation  effects  can  be  observed  in  which  the  net 
gaui  initially  overshoots  the  cavity  threshold  and  subsequently  decays  to  the 
steady  state  gain  condition  a{t)  = Under  conditions  cf  strong  pumping 

the  overshoot  can  be  significant,  resulting  in  an  initially  enhanced  super- 
radiative  extraction,  since  the  stimulated  emission  from  the  medium  occurs 
under  conditions  of  high  gain.  This  results  in  a sharpening  of  the  leading 
edge  of  the  optical  output  power,  and  is  observed  experimentally. 


3.3  Electron  Kinetics.  The  present  analysis  of  an  electrically  excited 
laser  mixture  features  a completely  coupled  description  of  the  molecular  and 
electron  kinetics.  Plasma  calculations  consist  of  numerical  solution  of  the 
Boltzmann  transport  equation  for  the  electron  energy  distribution  function, 
from  which  the  electron  transport  coefficients,  elastic  and  inelastic  collision 
rates,  and  electrical  power  partitioning  can  be  obtained  as  a function  of  gas 
mixture,  E/N,  and  degree  of  vibrational  excitation.  Superelastic  electron 
collisions  (i.  e.  , reverse  processes  in  which  electrons  gain  energy  from 
collisions  with  excited  molecular  species)  are  included  in  both  the  molecular 
and  plasma  kinetic  equations.  Basic  features  of  the  electron  kinetics  analysis 
are  summarized  in  Figure  33. 
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• Numerical  Solution  of  the  Boltzmann  Transport 
Equation;  Electron  Energy  Distribution  Fuiiction 

• Elastic  and  Inelastic  Processes  included:  Momentum 
transfer,  Vibrational,  Electronic  E-xcita+ion,  Ionization, 

. Attachment,  Reambination 

• Electron-electron  allisions  (important  for  fractional 
ionization  > 10  h 

• Superelastic  collisions  included  (necessary  for  analysis 
of  excited  gas  mixtures) 

• Electrical  parameters  (Mobility,  Drift  Velocity,  Average 
Energy,  Effective  Tern perstu re.  Plasma  Conductivity, 
Discharge  Current  Densih/,  etc. ) as  a function  of  E'M, 
Gas  Mixture 

• Fractional  Power  Partitioning  of  Electrial  Power  into 
all  elastic  and  inelastic  mechanisms 

• Forward  (and  reverse)  electrical  excitation  rates  <va> 
for  all  inelastic  collision  processes 


Figure  33.  Plasma  Kinetics  Analysis 
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In  general,  the  electron  energy  distribution  function  for  the  gas  mixtures 
and  relatively  low  E/N  values  typical  of  e -beam  sustained  electric  discharge 
lasers  is  highly  non-Maxwellian.  The  availability  of  extensive  experimental 
data  for  elastic  and  inelastic  electron  cross  sections  as  a function  of  energy 
for  many  gases  of  interest  often  makes  it  possible  to  implement  the  electron 
kinetics  analysis  by  solution  of  the  Boltzmann  equation.  The  electron  velocity 
distribution  function  f{v)  for  a multicomponent  gas  at  temperature  T in  an 
external  electrical  field  E satisfies  the  Boltzmann  transport  equation, 


[ 


ht 


+ V 


E (r,  t)  • ] f(v,  r,  t)  = -^ 


coll. 


(7) 


where  the  right-hand  side  collision  term  includes  the  effects  of  elastic  and 

inelastic  scattering  processes.  Discussion  of  the  Boltzmann  transport 

(41)  (42)  (43) 

equation  can  be  found  in  work  by  Holstein  , Allis  .Shkarofsky  et  al  , 

Frost  and  Phelps  ^'^^^,Carleton  and  Megill^^^\  and  Englehardt  et  (47)^ 

(41) 

T.  Holstein,  "Energy  Distribution  of  Electrons  in  High  Frequency  Gas 
Discharges",  Phys.  Rev.  70,  367  (1946). 

(42) 

' ' W.  P.  Allis,  "Motions  of  Ions  and  Electrons",  in  Handbuch  der  Physik. 

Vol.  XXI  (Springer-Verlag,  Berlin,  1956). 

I,  P.  Shkarofsky,  T.  W.  Johnston,  and  M,  P,  Bachynski,  The  Particle 
Kinetics  of  Plasmas.  Addison-Wesley  (1966). 

^44) 

L.  S.  Frost  and  A.  V.  Phelps,  "Rotational  Excitation  and  Momentum 
Transfer  Cross  Sections  for  Electrons  in  and  N2  from  Transport  . 
Coefficients",  Phys.  Rev.  127,  1621  (1962). 

(45) 

' ' N.  P.  Carleton  and  L.  R.  Megill,  "Electron  Energy  Distribution  in 

Slightly  Ionized  Air  under  the  Influence  of  Electric  and  Magnetic  Fields", 
Phys.  Rev.  126,  2089  (1962). 

A.  G.  Englehardt  and  A.  V.  Phelps,  "Elastic  and  Inelastic  Collision 
Cross  Sections  in  Hydrogen  and  Deuterium  from  Transport  Coefficients", 
Phys.  Rev.  131,  2115  (1963). 

G.  Englehardt,  A.  V.  Phelps,  and  C.  G.  Risk,  "Determination  of 
Momentum  Transfer  and  Inelastic  Cross  Sections  for  Electrons  in 
Nitrogen  Using  Tran  sport  Coefficients",  Phys.  Rev.  A135,  1566  (1  )64). 
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Various  approximations  and  techniques  for  numerical  solution  of  the  equation 
and  its  application  to  electric  discharge  gas  lasers  have  been  reported  by 
Nighan  ^,Rockwood^^*^^,Hancock  ^ Elliott  ^ ^^,Morgan  and 

Fisher^^^^jLowke  ^^^^^^,and  others. 


At  the  gas  densities  of  interest  for  laser  operation,  elastic  and  inelastic 
electron-molecule  collision  frequencies  are  so  high,  relative  to  molecular 
kinetic  processes,  that  the  electron  distribution  is  usually  adjusted 
instantaneously  on  the  time  scales  over  which  the  molecular  levels  are 
evolving.  Furthermore,  since  elastic  collision  frequencies  are  typically 
greater  than  inelastic,  the  electron  distribution  can  be  assumed  to  have  attained 
a quasi-steady  state  even  on  the  time  scales  of  inelastic  collision  processes. 

It  follows,  therefore,  that  the  first  tern,  of  (7)  containing  the  time  derivative 
can  be  discarded  if  the  electric  field  is  assumed  to  be  dc  or  slowly  varying. 

I or  a spatially  uniform  field,  it  is  similarly  possible  to  neglect  the  second 
term,  since  the  mean  free  path  (which  is  related  to  the  elastic  collision  cross 


(48) 

(49) 


W.  L.  Nighan,  "Electron  Energy  Distributions  and  Collision  Rates  in 
Electrically  Excited  N^,  CO,  and  CO^",  Phys.  Rev.  2A,  1989  (1970). 

S.  D.  Rockwood,  "Elastic  and  Inelastic  Cross  Sections  for  Electron-Hg 
Scattering  from  Hg  Transport  Data",  Phys.  Rev.  2348  (1973). 


J.  H.  Hancock,  R.  C.  Jones,  and  C.  B.  Mills,  "Numerical  Solution  of 
the  Boltzmann  Equation  for  Energy  Distributions  of  Electrons  with 
Inelastic  Scattering  on  Molecules",  Los  Alamos  Rept.  LA-4832,  May  1972. 


(51) 

C.  J.  Elliott,  O.  P.  Judd,  A.  M.  Lockett,  and  S.  D.  Rockwood, 

"Electron  Transport  Coefficients  and  Vibrational  Excitation  Rates  for 
Electrically  Excited  CO^  Gas  Lasers",  Los  Alamos  Informal  Rept. 
la-5562. MS,  April  1974. 

(52) 

W.  L.  Morgan  and  E.  R.  Fisher,  "Calculations  on  the  Electron  Energy 
Distribution  in  a Molecular  Laser  Plasma",  Research  Institute  for 
Eng.  Sciences,  Wayne  State  University,  Report  RIES  74-56,  August  1974. 

(53) 

J.  J.  Lowke,  A.  V,  Phelps,  and  B.  W.  Irwin,  "Predicted  Electron 
Transport  Coefficients  and  Operating  Characteristics  of  CO^-N  -He 
Laser  Mixtures",  J.  Appl.  Phys.  44.  4664  (1973). 
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section  and  molecular  density)  is  typically  small  compared  to  the  discharge 
dimensions.  The  Boltzmann  equation  (7)  then  becomes 


- (e/m)  E • 7 f(v)  = 6f/6t 

V 


coll 


(8) 


and  the  usual  approach  for  solving  this  equation  begins  with  an  expansion  of 
f(v)  into  Legendre  polynomials. 


f(v)  = f (v)  + (v/v)  • f (v)  + (3w  /v^  - 1):  f.  (v)  + . . . 
o 1 i. 


(9) 


(1  and  f^  represent  dyadic  quantities).  Usually,  f can  be  approximated  using 

only  the  first  two  terms.  Physically,  the  justification  of  a two-term 

approximation  is  related  to  assumptions  about  the  degree  of  anisotropy  to  be 

expected  for  the  velocity  distribution.  If  the  electric  field  is  small  enough 

that  the  directed  speed  of  the  electrons,  as  measured  by  their  drift  velocity, 

is  much  less  than  their  random  thermal  velocities,  the  small  first-order 

perturbation  v • f from  an  isotropic  distribution  f^  should  be  a good 

approximation.  For  high  values  of  E/N  characteristic  of  self-sustaining 

(54) 

discharges  and  for  certain  gases,  it  has  been  speculated  that  the  term 
may  be  comparable  to  the  f ^ term,  although  no  analysis  has  included  these 
higher  order  terms.  Retention  of  more  than  the  first  two  terms  requires 
knowledge  of  angular  (differential  scattering)  cross  section  data,  whicn  is 
not  generally  available.  Typical  experimental  data  gives  only  integrated 
cross  sections  as  a function  of  the  electron  energy.  Some  of  the  available 
data  in  the  literature  has  been  determined  indirectly  by  fitting  measured 
transport  coefficients  to  synthetic  cross  sections  using  a numerical  Boltzmann 
analysis  based  on  the  two-term  expansion  of  f(v).  Thus,  consistent  use  of 
such  data  requires  that  the  ele.tron  kinetics  calculations  should  be  based  on 


W.  F.  Bailey,  Paper  DD5  presented  at  the  Twenty-Fifth  Gaseous 
Electronics  Conf. , London,  Ontario,  Canada,  October  1972. 
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the  two-term  approximation.  For  values  of  E/N  typical  of  e -beam  sustained 
discharges,  this  approach  is  justified,  and  is  numerically  much  simpler  than 
would  be  the  case  if  angular  calculations  were  required. 


In  the  expansion  (9),  Ij('V’),  f^,  . •.  are  functions  only  of  the  magnitude 

V,  since  the  angular  dependence  occurs  in  the  spherical  harmonic  coefficients. 

It  can  easily  be  shown  that 


7 f = 
V 


_d_ 

dv 


(v  f , ) + V + 

dv 


V d 

3 dv 


(vf^)  • (3  vv  - 1)  -I-  . . . 


(10) 


where  f has  l.,'*en  separated  explicitly  into  terms  which  have  = 0,  1,  and  2 
symmetry,  respectively.  If  (10)  is  substituted  into  (8),  and  separated  into 
scalar  and  vector  parts  (which  corresponds  to  multiplying  by  1 or  v and 
integrating  over  the  spherical  solid  angles  Q^),  the  following  equations  ’•esult: 


eE  d , 2-?  . 6f 
7 • (v  f ) = o 

3m  V dv  6t 


coll 


eE 

m 


6t 


coll 


(11) 


(12) 


The  most  important  physical  mechanism  for  reducing  the  asymmetry  in  the 


distribution  function  is  electron-molecule  momentum  transfer  collisions.  The 

“*  (43) 

collisional  rate  of  change  of  fj^  can  be  approximated  by 


6t 


coll 


= -I  V ^ (v)  f 
X ™ ^ 


(13) 


X X 

where  = N^vQ^  (v)  is  the  momentum  transfer  collision  frequency  for 


molecules  X in  terms  of  the  population  density  N and  momentum  transfer 
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cross  section  Q . In  terms  of  the  electron  energy  u = mv  /2e  (in  units  of 
m o / » 

eV),  the  equation  for  f^  becomes,  after  manipulation  of  (11),  (12),  and 
substitution  of  (13), 


- (E  /3)  d/du 


-In 


u 





df  /du 
o 


]\P^= 


n/u 


6f 
c 

6t 


coll 


(14) 


Expressions  for  the  collision  term  on  the  right  hand  side  have  been  developed 
by  Allis  , Shkarofsky  ^ \ Frost  and  Phelps^^^  and  others  for  elastic 

momentum  transfer  and  inelastic  iiectron-neutral  scattering  processes.  In 
addition,  whenever  the  fractional  ionization  is  > lO"^,  it  is  often  necessary^^^^ 
to  include  the  effects  of  electron-electron  Coulomb  scattering.  Rockwood^"^*^^ 
has  shown  how  the  Boltzmann  equation  may  be  expressed  in  terms  of  a "flux 
divergent  description  of  the  electron  distribution  in  energy  space,  viz.. 


-dJ  /du  - dJ  /du  - dJ  /du  +Z[(u  + u ) f (u  + u ) N Q (u  4 u ) - 
^ cl  ee  Of  01  o Of  a a a 

a£^(u)N^*Q.^(u)]  = 0 


(15) 


where 


(u)> 

m 

J ,(u)=  -[u^  <(2m/M)NQ  > (f  + kT/e  df  /du)l 

Cl  mo 


4 

Je^(u)  = -(2TT/3)  ^ n^  In  A [P{u)  df^/du  4 Q(u)  f^(u)  ] 

U to 

P(u)  - 2 f dw  w^'^^  f (w)  4 [dw  f (w) 

J o ^ o 


(16) 

(17) 

(18) 

(19) 


(55) 


W,  H.  Long,  "Electron  Kinetics  in  the  KrF  Laser",  Appl.  Phys.  Lett 

31,  391  (1977). 
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(20) 


1/2 

Q(u)  = sjdww 
o 

In  the  definition  of  ^ = 300  e = 4.  8 x 10  esu  is  the  electronic  charge 

(cgs  units),  and  A = . , where  X and  r . are  the  Debye  lengths  and 

D min  D min 

classical  distance  of  closest  approach,  defined  by 


/ = (kT  /4TTn  ^)^^^ 
Dee 

r . = q^/(eu)  = q^/(l.  5kT  ) 

min  e 


(21) 

(22) 

(23) 


Q a.re  (forward  and  reverse)  cross  sections  for  the  inelastic  collision 

processes,  labeled  by  or,  and 

<NQ  (u)>  = (u) 

m ^ A m 

<(2m/M)  NQ  (u)>  =I  {2m /M„)  (u).  (24) 

ni  ^ A A m 

* 

and  are  population  densities  of  the  lower  and  upper  (excited)  states,  and 
u^  is  the  inelastic  energy  transferred  during  the  collision.  The  principle  of 
detailed  balance  provides  relations  between  the  cross  sections  for  forward 
and  reverse  processes: 


(u  + u ) Q (u  + u ) = uQ  (u) 

or  a a -a 


(u  - u ) Q (u  - u ) = uQ  (u) 
o -a  a or 


(25) 


The  fi:.st  two  terms  in  the  sum  over  or  in  (15)  correspond  to  collisions  in  which 
the  electrons  lose  energy,  while  the  last  two  terms  (which  represent  the 
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superelastic  processes)  correspond  to  collisions  in  which  excited  molecular 
species  transfer  energy  to  the  electrons.  Only  binary  electron-molecule 
collisions,  includmg,  for  example,  vibrational  excitation,  electronic 
excitation,  attachment,  recombination,  ionization,  etc.  will  be  considered. 

It  should  also  be  noted,  therefore,  that  the  collision  term  for  process  a 
does  not  always  contain  all  four  of  the  terms  indicated  in  the  general  expression 
given  in  (15).  For  ey  Tiple,  only  the  second  term  occurs  for  recombination 
or  attachment,  and  foi-  ionization  or  dissociation,  only  the  first  two  terms 
occur.  Effects  of  recoil  in  elastic  energy  transfer  pro  iesses  are  included  in  - 
the  second  term  of  (17),  where  the  molecular  distribution  is  assumed  to  be 
Maxwellian  at  temperature  T«  Effects  of  rotational  excitation  can  be  included 
in  the  same  way  ss  momentum  transfer  by  making  a suitable  continuum 
a p pr  oximation . 


The  normalizat’on  condition  for  f (u)  is  taken  to  be 

o 


f du  u^^^  f (u)  = 1 (26) 

J 0 

o 

1/2  2 

where  the  factor  u originates  from  the  density  of  states  v dv/du.  From 

the  solution  for  f (u),  all  of  the  forward  and  reverse  excitation  rates, 
o 

electrical  power  partitioning,  and  plasma  parameters  (drift  velocity,  mobility, 
average  and  characteristic  energies,  effective  temperature,  etc.)  can  be 
obtaLced.  The  important  quantities  required  for  the  coupled  laser  analysis 
are  the  secondary  electron  excitation  rates,  given  by 
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n (2e/m) 
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1/2 
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du  u f (u)  Q (u) 
o er 


(27) 


and  the  electrical  power  balance  (W/cm  ),  obtained  by  multipl', ing  (15)  by 
e(2e/rr.'^^^u  and  integrating: 
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efi  = e(2e/m)^^^ 


du  u < (2m /M)  NQ  > [f  + (kT/e)  df  /dul 
m o o'* 


o 


(28) 


where  |a  is  the  mobility.  The  physical  significance  of  (28)  is  that  the  total 
electrical  power  input  (on  the  LHS)  is  balanced  by  elastic  heating  from 
momentum  transfer  collisions,  and  by  power  dissipation  from  inelastic 
excitation  (and  de- excitation)  processes.  Note  that  there  is  no  contribution 
from  elastic  electron-electron  scattering:  this  mechanism  merely  redistributes 
electron  energy  and  tends  to  drive  the  distribution  toward  a Maxwellian  shape. 

Note  also  that,  for  processes  such  as  recombination  and  attachment,  in  which 
electrons  are  lost,  only  the  second  term  of  the  collision  expression  in  (15)  occurs, 
and  in  that  case,  the  corresponding  term  for  the  power  balance  (28)  must  be 
replaced  by 

00 

e (2e/m)^^^  N J du  u^f  (u)  Q^(u) 
a o Of 


Expressions  for  miscellaneous  plasma  parameters  are: 


p = - (1/3)  (2e/m^^^  [ du  (u/<NQ  (u)>)  df  /du 
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u = [ du  u^^^  f (u)  = (3/2)  kT  (29) 
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Numerical  solution  of  the  Boltzmann  equation  is  carried  out  by  reducing 
the  differential  equation  (15)  to  a finite  difference  equation,  defined  over  a 
uniform  grid  of  equally  spaced  energy  values.  This  results  in  a large  set 
(typically,  > 500)  of  coupled  algebraic  equations,  which  is  solved  by  an 
iterative  approach  that  reduces  the  system  to  a tridiagonal  set  of  equations. 

The  algorithm^^^^  which  has  been  implemented  is  extremely  fast,  and  has 
been  constructed  in  such  a way  that  execution  speed  is  optimized  when  a good 
ihitial  guess  is  available.  Thus,  the  present  approach  is  quite  suitable  to  a 
coupled  analysis  in  which  the  excited  state  population  densities  and  discharge 
voltage  change  as  a function  of  time,  since  the  previous  electron  distribution 
always  provides  a good  initial  guess  for  obtaining  an  updated  distribution. 

For  e.cample,  with  an  electron  grid  of  500  points,  15  inelastic  collision  processes 

and  electron- electron  scattering,  the  electron  distribution  converges  to  an 

-3 

accuracy  < 2 x 10  after  ~8  iterations  in  ~0.  7 CP  sec.  (Sample  output  is 
present  in  Figure  41,  42  and  43  in  the  following  section. ) 

3.4  External  Driving  Circuit.  Basic  features  of  the  equations  which 
describe  the  external  driving  circuit  are  illustrated  in  Figure  34,  A plane 
parallel  discharge  of  area  A,  with  anode  and  cathode  separated  by  a distance 
d,  is  driven  by  an  external  RLC  circuit  with  a capacitor  initially  charged 
to  high  voltage.  The  discharge  is  sustained  by  an  external  e -beam  with 
an  arbitrary  temporal  current  density.  Initially,  there  is  no  voltage  across 
the  discharge  because  of  the  external  circuit  inductance.  As  the  secondary 
electron  density  in  the  gas  changes,  the  plasma  conductivity  (and  hence 
the  discharge  impedance)  will  change.  Thus,  the  plasma  kinetics  analysis  is 
directly  coupled  to  the  circuit  equations  which  define  the  instantaneous 
voltage  (and  thus,  the  value  of  E/N)  across  the  discharge: 

^ ^ W,  B,  Lacina,  "Supersonic  Continuous  Wave  Carbon  Monoxide  Laser 

Development,  Volume  I:  CO  Laser  Kinetics  Code",  Northrop  Report 
NRTC-75-25R,  July  1975. 
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Plasma  conductivity;  o • ng(t)e^(t) 
Plasma  imoedanca:  Rj^lt)  * d/Aa(t) 

L li  + (r  + RJD)  i + 1 . 0 
dt  ^ C 


Figure  34.  External  Drivinc  Circuit 


L di/dt  + (r  + R (t))  i + q/C  = 0 
d 


(33) 


where  L = inductance,  r = external  resistance,  C = capacitance,  and  i = current. 
The  discharge  impedance  R^(t)  is  a function  of  time,  and  is  given  by 

r^(t)  = d/Aff{‘^),  (34) 

where 

ff(t)  = n (t)  e p (t)  (35) 

e 

and  p (t)  is  the  electron  mobility.  The  electron  kinetics  analysis  requires 
(E/N)(t)  = i(t)  R^(t)/(d/N). 

3.5  Typical  Analytical  Results.  Sample  output  from  a typical  computer 
run,  taken  from  a simulation  of  experimental  conditions  to  be  discussed  in  a 
subsequent  section,  will  be  presented  here  to  illustrate  tiie  features  of  the 
analysis. 

The  molecular  kinetic  processes  presently  included  for  the  KrF  system 
are  summarized  in  Figure  35,  which  is  a reproduction  of  the  input  data  deck 
which  defines  the  reaction  scheme.  Figure  38  is  a reproduction  of  output 
summarizing  details  and  comments  about  the  reaction  scheme;  this  table  is 
generated  during  the  synthesis  of  the  rate  subroutines  as  the  input  reaction 
scheme  is  processed,  scanned  for  errors,  and  translated  into  computer  coded 
equations.  Figure  37  presents  the  results  of  an  edit  of  the  reaction  scheme: 
the  occurence  of  each  species  encountered  in  each  reaction  is  summarized. 

As  Figure  35  and  36  indicate,  the  kinetic  processes  for  KrF  include  excitation 
of  excited  states  and  ions  of  Ar  and  Kr  by  secondary  and  high  energy  beam 
electrons,  fluorine  attachment,  formation  and  quenching  of  several  excimers, 
charge  transfer,  electron  and  negative  ion  (F  ) recombination  with  positive 
ions,  miscellaneous  replacement  reactions,  a variety  of  energy  transfer 
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SUMMARY  or  CARO  IhACES  TOR  INPUT  DATA  DECK 
(OATE!  10/21/771 


1 2 3 A 5 6 7 a 

123A567890123A567890123A567a90123A567a90123A567a90123A567890123A567890123A567890 


NORTHROP  KRF  LASER  KINETICS 
AR  ♦ E AR«  • E 


KR  ♦ E I.  KR*  » E 
AR  ♦ E ..  AR««  ♦ E 
KP  ♦ E ..  KR»«  • E 
aR«  ‘El*  aR«*  ♦ E 
KR«  ‘El*  kR«*  ♦ E 
AR*Ei*AR(*I  ‘E^E 
KR*Ef*KR(*)  «E«E 
AR»  •E.»AR(«)  ‘E^E 
KR«  ■•E(»kR(*)  ♦E«E 
aR««  ‘Ei^ARl^ltE^E 
kR«»  ‘El*  KR(.|  ♦ E ♦ E 
AR2<0  • E >•  AR*  * aR 
KR2<*»  ♦ E «•  KH«  ♦ KR 
F2  ♦ E ..  r « F- 


SCHAPERt  SCHElBNERt  BEITR.  AUS  PLASMA  PhTS.  9.  a5 
SCHAPERt  SCHEIBNERt  BEITR.  AUS  PLASMA  PHTS.  9,  A5 
0.1  X aR*  cross  SECTlONt  13. A EV  THRESH 
0.1  X aR*  cross  section.  11.5  EV  THRESH 
60.E-16  cm2  at  2 EV  ThRESHHOLO  ASSUMEO. 

60.E-16  CM2  AT  1.6  EV  THRESHHOLO  ASSUhED. 

RAPP.  ENGLANDER -COL DEN  JCP  A3.  1A6<>  (65) 

RAPP.  ENGlANOER-GOLDEN  JCP  A3.  1A6A  (65) 

7.E-16  CP2  AT  A.2  EV  THRESHHOLO  ASSUMED. 
lO.E-16  CM2  AT  A.O  EV  THRESHHOLO  ASSUMED. 

30  E*16  cm2  at  2. a EV  THRESHOLD. 

30  E-16-CM2  AT  3.1  EV  THRESHOLD. 

MEHR.  BIOnOI  (SYNTHETIC!  TE»»(-3/2)  LAO) 

OSKAM,  MITTELSTAOT  (SYNTHETIC!  TE»»(-3/2)  LAO) 

hoflano.  aerospace  CORP. 


aR2«  « C <•  AR  « AR  • E 

1.00  E-07 

kR2«  • E f KR  . kR  ♦ E 
1.00  E-07 

aR  • HE-  ..  AR(.)  ♦ he-  ♦ E 

A.3A  E-lO  BERGER-SElTZER:  STOP  * 1.7  MEVCM2/GM  j300  KEV) 

AR  ♦ HE-  ,.  aH»  ♦ HE- 

1.2a  E-ld  1.2a  E-18  (1/3.5)  X AR  . hE  AR(.)  ♦ E • HE 

KR  » HE-  ^ KR(.)  ♦ he-  ♦ E 


123A567890 1236567090 1236567890 123A 5678901236567890 123AS67B90123AS67890123A567890 
1 2 3 A 5 6 7 8 


Figure  35.  Summary  of  Card  Images  for  Input  Data  Deck 
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SUMMARY  OF  CARO  IMAGES  FDR  INPUT  DATA  DECK 
(OATe:  10/21/77) 


I 

r 

I 

I. 


CAHd  12  3 4 5 6 7 8 

NO.  l234S6  7e90l234S678901234S67e90|234S67a<»0 1234567890 1234567890 1234567890 1234567890 


a1  ....  8.04  E-18  BEBGER-SELT2Eflt  STOP  ■ 1.5  MEVCH2/GM  (300  KEV) 

42  ...,kR  * HE-  f.  KR»  . HE- 

43  ....  2.30  E-ie  2.30  E-18  (1/3.5)  A KR  . hE  <•  KR(.)  » E . HE 

44  ....AR(.)  . AR  . M I*  AR2(.)  • M 

45  ....  2.00  t-31  «.  F.  Lit’.  0.  C.  CON)*AY  JCP  62.  3070  (1975) 

46  ....AR(.)  • KR  . M ARKR(.)  . M 

47  ....  2.5E-31 

48  ....kR(*)  * KR  . M f.KR2(*)  ♦ M 

AV  ....  2.40  E-31  C.  J.  TRACY.  H.  J.  OSXAM  JCP  65.  3387  (1976) 

50  ....kR(*)  • AR  • M c ARKR(.)  » M 

51  ....  l.OE-31 

52  ....aR(«)  • KR  I*  AR  . KR(*) 

53  ....  3.00  E-11 

54  ....AR2(«I  * KR  f KR(*)  * AR  ‘ AR 

55  ....  7.5E-10  BOHME  ET  AL.  J.  CHEH.  PHYS.  52.  5094  (1970) 

56  ....ARKR(*)  • KR  >•  KR2(«)  . AR 

57  ....  3.2E-10  BohME  eT  aL.  J.  ChEM.  PHYS.  52.  5094  (1970) 

58  ....aR*  . AR  . M I*  AR2*  ♦ M 

59  ....  1.00  E-32  HILL.  GUTcHECk.  HUESTIS.  ET  AL.  SRl  REPORT.  1974. 

60  ....AR*  . KR  . M f.  ARKR*  . M 

61  ....  l.OE-32 

62  ....KR*  . KR  ♦ M <*  KR2*  . M 

63  ....  5.5E.32  HUGHES  LASL  ASPEN  9/76 

64  ,...KR*  • AR  . M I*  ARKR*  « M 

65  ....  I.OE-32 

66  ....ARkR*  ‘ KR  f.  KR2*  ♦ AR 

67  ....  I. 00  e-10  SRI  REPORT  NO.  HP  76-99.  DEC..  1976 

68  ....aR2<«)  ♦ F-  «•  ARF*  . AR 

69  ....  5.00  E-07 

70  ....AR2(«>  ‘ r-  » AR2F* 

71  ....  5.00  E-07 

72  ....AR(«>  ‘ F-  «•  ARF* 

73  ....  l.OE— 06 

74  ....kR2(‘I  ♦ F-  f.  KRF*  . KR 

75  ....  5.00  E-07 

76  ....kR2(*I  . F-  ..  KR2F* 

77  ....  5.00  E-OI 

78  ....KR(.)  ♦ F-  f.  KRF* 

79  ....  I.OE-06 

80  ....ARkR(*I  . F-  <•  KRF*  . AR 


1234567890 1234567890 1234567890 1234567890 123456789(1 1234567890 1234567890 1234567890 
1 2 3 4 5 6 7 8 


Figure  35  (continued) 
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SUMMARY  or  CARO  iMACCS  FOR  INPUT  QATA  DECK 
(OATEi  10/21/771 


CAHO  12  3 4 5 6 7 6 

NO.  12345678901234567890123456789012345678901234567890123456789012345678901234567890 


81  ....  5.00  E-07 

82  ....ARKR<*)  ♦ F-  I*  ARKRF* 

83  ....  5.00  t-07 

84  ....kR«  » F2  » KRF«  • F 

85  ....  7.20  E-lO  VELAZCOf  KOLTSt  SETSERt  JCP  65t  3469  (1976) 

86  ....AR«  ♦ F2  c ARF*  * F 

87  ....  7.50  E-lO  VELAZCOf  KOl.TSt  5ET5ER.  JCP  65t  3469  (19761 

88  ....AR«  * Kr  ..  AR  * KR« 

89  ....  *-.26  E-12  PlPERf  SETSERt  CLYNEt  JCP  63t  4018  «1975) 

90  ....KRF*  ♦ AR  * M ft  ARKRF*  » M 

91  ....  6.00  E-32  AVCO 

92  ,...KRr«  • KR  » M ft  KR2F«  t M 

93  ....  5.00  E-31  AvCO 

94  ....ARF*  » AR  t M ft  AR2F*  » M 

95  ....  4.00  E-31  AVCO 

96  ....ARf*  • KR  t M ft  ARKRF*  ♦ M 

97  ....  1.00  E-31 

9B  ....AR2*  • KR  ft  AR  ♦ AR  t KR* 

99  ....  8.00  E-11  Zamir  (priv.  commun.  to  srii 

100  ....AR2*  * AR2*  ft  AR2<tl  . AR  . aR  * E 

101  ....  3.0E-10 

102  ....ARf*  • KR  ft  KRF*  ♦ AR 

103  ....  1.50  E-10  SRI  REPORT  NO.  MP  76-99.  DEC..  1976 

104  ....AR2*  * F ft  ARF*  . AR 

105  ....  3.0E-10 

106  ....KR?*  * F ft  KRF*  » KR 

lOT  ....  3.0E-10 

108  ....AR2*  ‘ F2  ft  aR2F*  . F 

109  ....  2.60  E-lO  SRI  REPORT  NO.  MP  76-99.  OEC..  1976 

no  ....aRKR*  • F2  f KRF*  . AR  » F 

111  ....  6.00  E-10  SRI  Report  no.  mp  7b-99.  oec..  1976 

112  ....ARKR*  • F2  -•  ARKRF*  t F 

113  ....  3.0E-10 

114  ....KR2*  ♦ F2  ft  KR2F*  » F 

115  ....  3.00  E-10  SRI  Report  no.  mp  76-99.  oec..  1976 

116  ....AR2F*  » F2  >•  AR  . AR  » F . F2 

117  ....  1.00  E-09  SRI  REPORT  NO.  MP  76-99.  OEC..  1976 

118  ....aRkRF*  . F2  ft  AR  ♦ KR  t F . F2 

119  ....  1.00  E-09  SRI  report  NO-  MP  76-99.  OEC..  1976 

':20  ....KR2F*  • F2  ft  KR  . KR  ♦ F t F2 


12345678901234567890123456789012345678901234567890123456789012345678901334567890 
1 2 3 4 5 6 7 8 


Figure  35  (continued) 
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SUMMARY  OF  CARO  IMAGES  FDR  INPUT  DATA  DECK 
ICATeI  10/21/77) 


card  12  3 4 5 6 7 8 

MO.  1234567890 1234S67e90 1234567890 123<>S678<)0 1234567890 1234567890 123456  7890 1234567890 


121 

....  1.00  E-09 

SRI 

repdht 

ND. 

MP 

76-99* 

DEC.* 

1976 

122 

....ARKRF*  . KR  I.  KR2F»  • AR 

123 

....  1.00  E-10 

SRI 

hepdrt 

NO. 

HP 

76-99* 

DEC.* 

1976 

124 

....KRf*  * F2  KR  * F . F2 

125 

....  1.00  E-09 

SRI 

repdrt 

nd. 

HP 

76-99. 

DEC., 

1976 

126 

....ARF*  * F2  - AR  * F . F2 

127 

....  1.00  E-09 

SRI 

repdrt 

nd. 

HP 

76-99* 

DEC.* 

1976 

128 

....AR2F4  ♦ KR  ..  KrF*  . AR  » 

AR 

129 

....  1 . OE- 1 0 

130 

....AR4*  • M ar»  ♦ m 

131 

....  1.00  E-10 

132 

....KR4«  • m ,.  kr»  • m 

133 

....  1.00  E-10 

134 

....ARF*  <•  AR  . F 

13S 

....  3.30  E 0/ 

SRI 

repdht 

nd. 

HP 

76-99* 

DEC., 

1976 

136 

....AR2«  1*  AR  . AR 

137 

....  3.80  E 06 

SRI 

repdrt 

ND, 

HP 

76-99* 

DEC.* 

1976 

138 

....ARkR4  <•  AR  . Kr 

139 

....  3 . OE  *06 

140 

....KR2*  f*  KR  . KR 

141 

***.  3*30  E 06 

SRI 

repdrt 

NO. 

MP 

76-99* 

DEC.* 

1976 

142 

....AR2F«  (•  AR  . ah  . F 

143 

....  2.00  E 08 

SRI 

repdrt 

nd. 

HP 

76-99. 

dec.. 

)976 

144 

....ARkRF*  ..  AR  . kr  ♦ F 

145 

....  5.00  E 07 

5R1 

report 

nd. 

HP 

76-99* 

DEC.* 

1976 

146 

....kR2F»  >•  kr  * kr  * F 

147 

....  6.70  E 07 

SRI 

repdrt 

nd. 

HP 

76-99* 

DEC.  ^ 

1976 

1^*8 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 


R,  BURNHAM*  S.  5.  SE ARLES  (SUBMITTED  TD  JCP) 


KRF»  »•  KR  • F •*  HNU 
1.10  b 08 
KRF*  ♦ RAO  I*  KR  ♦ F • RAD 

2.00  E-16  laser  TRANSITIONS  STIMULATED  EMISSION  X-5ECT1DN 

F2  ♦ RAO  ..  F ♦ F 

1.50  E-20 

F-  ♦ RAO  f*  F ‘ E 

5.40  E-18  a.  MAnDL*  PHY5  REV  A3*  251  (1971) 

KR2E*  • RAD  >•  AR*  • KR  ♦ F 

1.00  E-99  cross  SECTION  UNKNOWN 

AR2(*  I ♦ RAD  -•  AR  ♦ AR(*) 

1.50  E-17  STEVENS  (PARK  CITY  CONFERENCE) 

KRa(*)  • rad  KR  • KR(*) 


12345678901 234567890 123A567890 12345678901 2345678901 234567890 1234567890 123A567890 
1 2 3 4 5 6 7 8 


Figure  35  (continued) 
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Figure  36  (continued) 
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Figure  38.  Summary  of  Card  Images  for  Input  Data  Deck 
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Figure  38  (continued) 
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collisions  and  quenching  processes,  spontaneous  radiative  decay,  stimulated 
emission,  and  absorption  from  and  various  excited  transient  species  in  the 
active  medium.  Figure  38  illustrates  the  input  data  deck  for  execution  of  the 
analysis  for  conditions  typical  of  the  scale-up  device  experiments.  Figure  39 
summarizes  the  input  reaction  scheme  after  rate  constants  (if  any)  have  been 
modified  (none  for  the  present  example). 

Figures  40-57  illustrate  typical  output  generated  during  execution  of  the 
analysis.  The  formating  is  explicit  enough  that  the  material  presented  in 
these  figures  should  be  self-explanatory.  In  order  to  provide  some  quantitative 
assessment  of  the  sensitivity  of  the  various  reactions  included  in  the  kinetic 
scheme,  the  computer  code  generates,  according  to  output  request,  tables 
such  as  those  presented  in  Figure  45.  These  tables  summarize  the  numerical 
value  of  the  contribution  each  reaction  makes  (at  any  given  time)  to  each 
species  in  the  system,  expressed  for  each  species  as  a percentage  of  the 
maximum  ’•ate  occuring  for  that  species.  Those  reactions  which  never 
contribute  more  than  a specified  percentage  are  flagged  with  an  asterisk,  and 
are  summarized  at  the  conclusion  of  the  analysis.  Thus,  after  a variety  of 
calculations  have  been  made  for  different  experimental  conditions  and  parameters, 
it  is  possible  to  determine  whether  any  given  reaction  may  be  considered  to 
be  relatively  unimportant  and  deleted  from  the  kinetic  scheme. 

Figures  46-57  represent  a portion  of  the  computer  plots  generated  at  the 
terminatioii  of  the  analysis.  Included  are  temporal  plots  of  e -beam  current 
density,  radiation  intensity,  intracavity  gain  and  absorption  coefficients, 
electrical  inp  it  and  optical  extraction  power  densities,  discharge  current 
density  and  voltage,  plasma  conductivity  and  impedance,  output  energy  and 
efficiency,  and  (requested)  population  densities  of  excited  species. 


Due  to  multiple  scattering  effects,  energy  deposition  in  a high  pressure  gas  is 
somewhat  greater  than  that  given  by  the  Berger-Seltzer  cross  sections.  The 
effective  cross  sections  used  in  the  analysis  are  taken  to  be  a factor  of  2.0  higher 
(cf.  DEPOSIT  = 2.0  in  Figure  39). 
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Figure  39.  Sumrrnry  of  Updaled  Rales  for  Input  Reaction  Scheme 

of  Synthetic  Kinetics  Code 
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Figure  39  (continued) 
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Figure  40.  Population  Densities  and  Rates  of  Change 
Electrical  and  Optical  Parameters 
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Figure  41.  Plasma  Parameters 
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Figure  44.  Population  Densities  and  Rates  of  Change 
Electrical  and  Optical  Parameters 
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TINE  I > S.OOOE-ttT  SEC 

percentage  contribution  or  reaction  k to  onii»/ot«  expressed  iroR  each  speciesi 
AS  A percentage  or  the  maxinuh  rate  occijring  roR  all  reactions  included 


Figure  45.  Percenlage  Contrlbulion  of  Reaction  K to  DN(I)/DT 


Figure  45  (continued) 


THIS  reaction  contributes  less  Than  s.  « to  all  shecIES  throughout  the  entire  calculation  so 
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AS  A percentage  of  the  HAXIHUH  rate  OCCURING  for  all  reactions  ITTCLUDEO 
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TIHE  T • S,000E-07  SEC  . _ , 

percentage  contribution  of  reaction  k to  dnii)/dT>  expressed  troP  Each  species) 
AS  A percentage  of  the  NAAIHUH  RATE  OCCuRING  FOR  ALL  REACTIONS  INCLUOEO 
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• THIS  reaction  contributes  less  than  5.  % TO  ALL  SPECIES  THROUGHOUT  THE  ENTIRE  CALCULATION  SO  FAR 


Figure  45  (continued) 
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percentage  contribution  or 
AS  A percentage  or  THE 


TINE  T ■ S,000E-07  SEC 

REACTI'R*  K to  DN(l)/DTt  EXPRESSED  (roR  EACH  SPECIES) 
MAXIMUM  RATE  OCCURING  roR  ALL  REACTIONS  INCLUDEO 
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THIS  reaction  contributes  less  ThAK  S.  ft  TO  ALL  SPECIES  TI^OUGHOUT  THE  ENTIRE  CALCULATION  SO  EAR 


Figure  45  (continued) 


-BEAM  CUBHENi  OENSIIV  IAH|>/CH?»  AS  A FUNCIION  OF  TIHE 


Figure  46.  E-Beam  Current  Density  as  a Function  of  Time 


Figure  47.  I nlracavlly  Radiation  I nlenslly 


Figure  49.  Laser  Gain  and  Medium  Absorption  Coefficients 
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Iqure  52.  Optical ly  Extracted  Power  Density  and  Efficiency 
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Figure  54.  Plasma  Conducllvlty  and  Discharge  Impedance 
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4.  EXPERIMENTAL  RESULTS  AND  THEORETICAL  CORRELATION 


This  section  contains  a detailed  presentation  of  experimental  results  and 
a discussion  of  device  characteristics  which  limited  the  accessible  range  of 
experimental  parameters.  Comparisons  of  these  results  with  the  theoretical 
model  has  been  undertaken  and  the  correlation  has  been  generally  good  in  all 
essential  details.  Results  of  a sample  comparison  are  presented.  In 
addition,  an  analysis  is  presented  which  attempts  to  explain  limitations  in 
the  experimental  results  based  on  nonuniformity  in  electron  beam  energy 
deposition.  Finally,  a theoretical  extrapolation  is  presented  assuming 
optimistic  device  characteristics. 

The  experimental  results  fall  basically  into  two  classes  dictated  by  the 
electron  beam  deposition  uniformity.  The  first  class  deals  with  laser 
extraction  over  small  volumes  with  the  discharge  anode  placed  3-5  cm  from 
the  ground  plane  screen.  In  this  case,  total  energy  output  was  sacrificed  in 
an  attempt  to  determine  discharge  stability  limits  and  enhancement  with 
reasonable  beam  deposition  uniformity  over  short  distances.  In  the  second 
case,  the  discharge  anode-cathode  separation  was  set  at  its  extreme  value  of 
10  cm  with  the  intent  of  extracting  maximum  energy  at  the  expense  of  discharge 
enhancement.  In  both  instances  of  small  and  large  anode-cathode  separations, 
for  e -beam  pumping  alone,  the  total  energy  output  was  found  to  be  propor- 
tional to  electron  beam  current  density,  whereas  energy  enhancement  v'Uh 
the  addition  of  a discharge  was  found  to  be  roughly  inversely  proportional  to 
this  parameter.  Here,  enhancement  is  defined  as  the  ratio  of  laser  energy 
output  with  electron  beam  plus  discharge  excitation  to  laser  energy  output 
with  electron  beam  excitation  only.  Under  optimum  operating  conditions 
with  a 900  ns  electron  gun  pulse  length,  the  maximum  extracted  specific 

energy  was  22  J/liter  for  a 3 liter  volume  with  an  enhancement  of  1.22  and 

2 

beam  current  of  18  A/cm  . Total  extracted  energy  at  this  current  level  was 
90  J over  10  liters.  For  a reduced  beam  current  density  of  10  A/cm^, 
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specific  extraction  was  likewise  reduced  to  8 J/liter  with,  however,  a 
corresponding  increase  in  enhancement  to  1.8.  The  relatively  low  enhance- 
ment factors  are  shown  below  to  be  a direct  result  of  discharge  instability 
induced  by  nonuniform  beam  deposition  in  the  direction  of  propagation. 

4.  1 Experimental  Results.  The  experimental  results  are  presented 
sequentially  in  four  basic  groups:  a qualitative  description  of  beam 
uniformity,  resulting  laser  output  patterns  and  discharge  profiles;  an 
analysis  of  typical  waveforms  for  the  KrF  and  XeF  systems;  a detailed 
outline  of  parametric  studies  to  determine  nominal  operating  conditions; 
and  a summary  of  results  under  optimum  conditions. 

The  question  of  electron  beam  deposition  was  fundamental  to  the 
experimental  results  and  has  been  detailed  in  Section  2.  1.  However,  a 
graphic  representation  would  be  helpful  here  and  is  shown  in  Figure  58. 

The  photograph  is  an  end  view  of  the  device  looking  down  the  optical  axis 
and  shows  simply  total  fluorescence  generated  by  the  electron  beam 
propagating  from  left  to  right  into  2 atm  of  Argon.  The  anode  was  drawn 
back  15  cm  and  can  be  distinguished  by  the  slight  line  reflection  at  the 
right  side  of  the  photograph.  The  gun  cathode  in  this  case  was  of  the  high 
current  density  variety  with  a 20  cm  width.  The  extent  of  maximum 
beam  fluorescence  is  about  17  cm  in  the  transverse  direction  and  about 
10  cm  in  the  direction  of  propagation. 

A typical  laser  output  burn  pattern  on  polaroid  film  is  shown  in  Figure 
59  for  electron  beam  excitation  only,  and  correspondence  with  the  previous 
measurement  of  energy  deposition  uniformity  is  apparent.  Energy  density 
decreases  from  the  center  outward  with  the  dark,  oval  center  section  containing 
the  highest  density  and  the  surrounding  light  halo  about  1/3  as  much.  At  the 
border  of  this  second  area  is  a third,  faintly  discernible  in  the  photograph. 
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(12  liters' 


10  cm 
(14  liters) 


Figure  59.  Laser  Output  Pattern 


which  contributes  very  little  to  total  energy.  The  rectangular  border 
indicates  the  position  of  the  optic  holder  which  extends  into  the  laser  plenum 
providing  for  an  optical  cavity  length  of  1.  5 M,  A support  web,  shown  at  the 
center  of  the  optic  holder  and  burn  pattern,  obscures  about  0.2  liters  of  gain 
volume.  The  burn  pattern  corresponds  approximately  to  a 10  liter  gain 
volume,  and  the  12  and  14  liter  designations  below  their  respective  lineal 
dimensions  indicate  the  available  extraction  volume. 

A"  mentioned  in  Section  2.  1 wherein  the  foil  support  geometry  was 
detailed,  major  support  web  obscuration  was  appreciable.  Under  discharge 
current  avalanche  conditions,  the  discharge  was  found  to  collapse  into 
multiple,  fairly  diffuse  channels  coincident  with  the  inter-web  open  spaces 
with  little  or  no  current  drawn  at  the  web  areas  themselves.  This  phenomenon 
is  shown  in  the  photographs  of  Figure  60.  As  indicated  in  Figure  60a,  the 
camera  was  tilted  to  observe  the  discharge  down  the  full  extent  of  the  anode. 
Figure  60b  shows  the  discharge  collapse  for  applied  fields  just  below  full 
breakdown  values.  The  central  image  is  the  one  of  interest;  the  images  to 
either  side  of  this  are  reflections  from  the  optic  extension.  The  discharge 
anode  is  to  the  left  and  gun  foil  to  the  right.  For  much  higher  applied  fields, 
the  discharge  breaks  down  fully  as  shown  in  Figure  60c,  and  the  most  intense 
luminosity  is  associated  with  the  region  of  initially  highest  beam  deposition. 
Since  no  sign  of  discharge  activity  larger  than  these  dimensions  was  observed 
on  the  anode  surface,  it  was  difficult  to  estimate  the  distribution  of  discharge 
current  in  the  more  expanded,  stable  regime.  Most  likely,  discharge  current 
followed  the  electron  beam  deposition  profile  in  which  case  it  would  closely 
match  the  beam  fluorescence  or  output  burn  patterns. 

The  areas  of  current  concentration  caused  blemishes  on  the  anode  as 
shown  in  Figure  61.  The  dark  spots,  roughly  2.  5 cm  wide  by  6 cm  long, 
coincided  exactly  with  the  inter -web  open  areas  and  the  luminous  channels 
of  Figure  60. 
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Figure  61.  Arc  Blerishes  on  Anode 


The  action  of  the  discharge  channels  on  the  anode  surface  and  the  bright 
arcs  in  the  photographs  were  a result  of  long  term  oscillatory  currents,  since 
after  electron  beam  turnoff  the  discharge  was  found  to  ring  for  times  on  the 
order  of  20  |x  s as  shown  in  Figure  62.  The  period  of  interest  was  in  the  first 
microsecond,  a small  fraction  of  the  total  waveform.  Therefore,  the  effects 
of  discharge  collapse  pr  ented  above  correlate  strongly  with  the  period  after 
termination  of  the  electron  gun  pulse  when  discharge  collapse  would  be 
expected  in  any  case  without  the  stabilizing  influence  of  the  electron  beam. 

These  arc  effects  were  not  observable  during  the  gun  pulse  owing  to  their 
small  contribution  to  the  temporal  averaging  measurements  employed  here. 
However,  the  long  term  current  spatial  behavior  is  useful  since  it  is  most 
likely  indicative  oi  that  occurring  during  the  gun  pulse  if  discharge  collapse 
occurs.  It  is  reasonable  to  assume  that  the  avalanche  distribution  resulting 
from  discharge  collapse  during  the  gun  pulse  would  carry  over  to  later  times 
as  well.  At  no  time  was  a distribution  other  than  that  outlined  above  observed, 
suggesting  that  the  arc  geometry  during  electron  beam  operation  was  as 
indicated  by  the  long  term  measurements. 

Figure  62  also  gives  an  indication  of  the  inherent  circuit  inductance  in  the 
case  of  discharge  collapse.  From  the  usual  relationship  relating  capacitance 
and  inductance  to  ringing  frequency,  the  inductance  is  determined  to  be  about 
225  nH.  When  compared  to  the  crowbar  inductance  of  175  nH  it  is  clear  that 
the  arc  channels  are  capable  of  competing  favorably  with  the  crowbar  circuit 
for  energy.  Ind'ed,  it  was  found  that  the  crowbar  did  not  alter  the  waveforms 
obtained  without  it. 

Typical  V/aveforms  . The  discharge  and  laser  output  were  characterized 
with  the  two  basic  circuit  types  reviewed  previously.  These  consisted  of  a 
simple  capacitive  discharge  and  a pulse  charged,  two  stage  capacitive  discharge. 
In  the  former  case,  it  was  found  that  significant  energy  v/as  delivered  only  in 
the  latter  half  of  the  excitation  pulse,  limited  by  discharge  cing.  With  the 
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Figure  62  Long  Term  Discharge  Voltage  and  Current 


pulse  charged  discharge,  it  was  observed  that  energy  was  delivered 
significantly  earlier  in  time,  but  that  discharge  arcing  also  occurred  earlier, 
with  the  net  result  that  the  two  circuits  were  equivalent  as  far  as  total  energy 
loading  was  concerned.  By  varying  storage  capacitance,  attempts  were  made 
to  reduce  discharge  voltage  late  in  the  excitation  pulse  in  order  to  prevent 
arcing  at  these  times.  It  was  observed  that  total  energy  input  did  not  vary 
significantly  for  capacitive  storage  in  the  range  1.  33-2.33  pF.  Above  these 
values,  arcing  was  not  observed  to  occur  sooner,  but  the  excess  stored 
energy,  channeled  into  several  columns,  would  often  lead  to  gun  foil  rupture. 
The  relevant  discharge  and  laser  output  waveforms  are  presented  in  the 
following  eight  figures  for  both  the  KrF  and  XeF  systems. 

Figure  63  shows  typical  traces  for  the  capacitive  discharge  with  an 
anode-cathode  separation  of  10  cm  and  central  beam  current  density  of  10  A/cm 
A nominal  gas  mixture  of  2 Torr  F^,  100  Torr  Kr  and  1,  7 atm  Ar  was  used. 
Voltage  values  along  the  top  of  each  column  indicate  initial  capacitor  bank 
voltage  and  are  included  for  labeling  purposes.  Relative  energy  output  is 
indicated  along  the  bottom  of  each  column.  Laser  output  with  beam  only 
excitation  is  shown  in  the  far  left  column  and  relative  output  in  this  case  is 
67%  of  the  ma-rimurn  obtained  with  discharge  excitation.  The  upper  trace 
displays  gun  voltage,  which  closely  approximates  the  current  waveform,  and 
the  lower  trace  indicates  laser  intensity  measured  with  a photodiode.  The 
second  column,  indicated  by  an  initial  charge  voltage  of  30  KV,  shows 
discharge  voltage  in  the  fairly  erratic  upper  trace  and  discharge  current  in 
the  smooth,  monotonically  increasing;  lower  trace.  The  discharge  was 
triggered  on  at  the  same  time  as  the  electron  gun;  therefore,  the  initial 
discharge  voltage  greatly  exceeds  the  nominal  trace  values  since  the  discharge 
impedance  is  high  at  early  times  before  appreciable  gun  current  is  established. 
The  possibility  that  high  discharge  field  strengths  prior  to  gun  turn  on  could 
set  up  incipient  arc  channels  which  would  lead  to  arcing  during  the  pulse  was 
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Figure  63.  Capacitive  Discharge 


tested  by  delaying  the  discharge  with  respect  to  gun  turn  on.  Identical  results 
were  obtained  for  all  discharge  delays  indicating  that  this  possibility  was  not 
realised.  After  full  beam  current  is  established  the  discharge  impedance 
drops  and  the  applied  discharge  voltage  is  likewise  lowered  at  which  point 
discharge  current  begins  to  flow.  As  the  rate  of  rise  of  discharge  current 
decreases,  the  inductive  voltage  drop  lessens  and  discharge  voltage  begins  to 
rise.  However,  voltage  begins  to  fall  off  roughly  500  ns  into  the  pulse  as  the 
capacitor  bank  is  depleted  until  finally  at  about  900  ns  the  electron  gun  turns 
off  and  the  discharge  impedance  rises  once  more.  This  results  in  a second 
sharp  peak  in  discharge  voltage  due  to  residual  energy  in  the  capacitor  bank. 
After  this  point,  voltage  dtcps  drastically  and  current  continues  to  rise 
indicating  the  onset  of  arcing.  As  the  discharge  waveforms  show,  discharge 
energy  is  pumped  in  primarily  late  in  the  pulse  and  this  is  reflected  in  the 
late  enhancement  of  laser  output.  These  trends  are  apparent  in  the  other 
traces  for  40  and  50  KV  initial  charge  voltages  as  well.  At  higher  charge 
voltages,  however,  discharge  voltage  drop  is  experienced  earlier  in  the  pulse, 
whereas  current  continues  to  rise  indicating  that  arcing  occurs  earlier  also. 

The  fact  that  the  laser  pulse  is  also  terminated  at  earlier  times  supports 
this  observation.  Unfortunately,  the  onset  of  arcing  cannot  be  more  firmly 
established.  The  usual  sign  in  higher  impedance  discharges  of  an  inflection 
in  discharge  current  is  lacking  here  due  to  the  low  discharge  impedance  which 
does  not  change  appreciably  under  discharge  collapse.  Attempts  at  observation 
of  increases  in  discharge  fluorescence  with  a photodiode  at  the  time  of  arcing 
failed,  since  fluorescence  levels  were  not  noticeably  different  before  and 
after  the  indicated  time  of  discharge  collapse  during  the  pulse. 

Typical  traces  obtained  with  the  pulse  charged,  two  stage  discharge  are 
shown  in  Figure  64;  and  the  various  parametric  values  were  identical  to  the 
pre^dous  case.  With  reference  to  the  20  KV  initial  charge  case,  the 
discontinuous  upper  trace  indicates  discharge  voltage  and  the  lower,  almost 
flat  curve,  indicates  capacitor  bank  current.  The  initial  steep  rise  in  discharge 
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Figure  64.  Pulse  Charged-Two  Stage  Discharge 


voltage  results  from  charging  of  C^,  originally  at  ground  potential.  Coincident 
with  this  is  a slight  rise  and  fall  in  current  from  more  easily  discerned 
in  the  higher  voltage  traces.  At  600  ns  into  the  charge  cycle,  the  electron 
beam  is  switched  on  and  the  discharge  impedance  collapses  resulting  in  a 
corresponding  drop  in  discharge  voltage.  After  this  point,  the  current  traces 
do  not  reflect  true  discharge  current  and  should  be  ignored.  Subsequent 
to  beam  turn  on,  the  discharge  voltage  behavior  is  identical  to  the  previous 
case.  It  is  observed  that  laser  intensity  is  enhanced  at  earlier  times  compared 
to  the  simple  capacitive  case,  but  that  termination  of  the  pulse  occurs  earlier 
at  the  higher  voltages.  This  effect  is  most  likely  due  to  the  fact  that  arcing 
field  strengths  are  approached  more  quickly. 

As  shown  in  Figures  63  and  64,  the  enhancement  ratio  for  the  simple 
capacitive  discharge  was  about  1.5  whereas  that  for  the  pulse  charged,  two 
stage  discharge  was  1.25.  The  difference  is  within  the  experimental 
variability  over  a period  of  time  for  a number  of  such  experimental  runs. 

These  previous  results  were  obtained  for  a non-optimized  system  and  are 
meant  to  represent  typical  discharge  behavior.  In  both  cases  cited  above, 
maximum  extracted  energy  density  was  7 J/liter  at  about  10%  efficiency. 

One  of  the  more  important  experimental  parameters  was  partial 

pressure,  and  results  of  variation  in  this  quantity  over  the  range  1.  5-3.0 

Torr  T.re  shown  in  Figures  65,  66  and  67.  In  this  series,  anode-cathode 

spacing  was  reduced  to  5 cm,  capacitive  storage  reduced  to  1.33  pF,  and 

/ 2 

beam  current  density  increased  to  15  A/cm  . For  the  1.5  Torr  case  of 
Figure  65,  almost  no  enhancement  was  obtained,  whereas  at  the  higher 
partial  pressures,  enhancement  was  about  1.23.  Also  for  increasing  partial 
pressure,  discharge  impedance  is  observed  to  increase  as  a result  o^ 
increasing  F^  attachment.  Highest  energy  extraction  was  obtained  for  2 Torr 

F^  which  was  considered  the  nominal  value,  and  output  for  the  1.5  Torr  and 
3 Torr  cases  was  down  by  10  and  20%  respectively  from  this  maximum. 
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Figure  65.  Typical  Waveforms -1. 5 Torr  F2 


Figure  66.  Typical  Waveforms-^  0 Torr  F2 


Experiments  were  performed  for  3 atm  total  pressure  and  reduced 
anode-cathode  spacing  in  order  to  ensure  beam  deposition  uniformity.  Typical 
results  are  shown  in  Figure  b8.  Characteristic  of  this  regime  is  the  spiked 
Isser  r'utput  waveform  which  was  exhibited  for  a range  of  gas  mixtures. 

XeF  Extraction  Resiolts.  Lasing  with  the  XeF  system  was  attempted  for 

(57) 

gas  mixtures  of  current  interest  involving  NF^/Xe/Ne  in  the  ratios 
1.67/5,0/1267  Torr  and  2,5/7.5/2280  Torr  and  the  results  are  shown  in 
Figures  69  and  70,  respectively.  Electron  beam  current  density  was  15  A/cm 
for  about  a 600  ns  pulse.  For  the  data  of  Figure  69,  a capacitive  discharge 
was  used  resulting  in  a specific  extraction  of  1.85  J/liter  at  about  1.0% 
efficiency  and  an  enhancement  of  2.2. 


Operation  at  higher  pressure,  as  indicated  in  Figure  70,  resulted  in  an 

increase  of  specific  extraction  to  2.5  J/liter  at  an  enhancement  of  1.5.  In 

this  latter  case,  the  pulse  charged,  two  stage  discharge  was  used  at  an  anode- 

cathode  separation  of  3 cm.  These  results  compare  favorably  with  those  rf 

others,  where  maximum  extraction  energies  of  2,8  J/liter  at  1,8%  efficiency 

(57) 

have  been  obtained  for  1 )Ji  s pulse  lengths  . 


Parametric  Curves.  Device  parameterization  for  gas  mixture,  total 
pressure  i.nd  discharge  voltage  are  presented  in  the  following  series  of 
curves.  In  order  to  guarantee  reliable  electron  gun  operation  for  the  complete 
set  of  data,  current  density  was  reduced  to  10-12  A/cm^  for  an  800  ns  pulse 
length.  Since  comparison  of  results  at  high  and  low  pressure  was  of  interest, 
the  discharge  anode-cathode  separation  was  reduced  to  6 cm  in  order  to 
provide  for  uniform  deposition  over  the  pressures  of  interest.  Discharge 
energy  was  supplied  by  the  simple  capacitive  circuit  with  = 2.  33  |iF;  and 


(57) 


L.  F.  Champagne  and  N.  W.  Harris,  Appl,  Phys.  Lett.  3^,  513 
(1977). 
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Figure  68.  Typical  Waveforms  at  High  Pressure 
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Figure  69.  Typical  Waveforms-XeF  at  Low  Pressure 


1 

1 


44 


A-K»3  cm,  V*l.  0 liter 
C|*i.  33 /iF,  C2*0. 7 ^fF 
NF3/Xe/Ne«k.  5/7. 5/2280  Torn 


2 KV/Div 
4. 5 KA/Div 


Laser  Intensity 
a5  J/liter 

Enhancement=l.5 


200  ns 


Figure  70.  Typical  Waveforms-XeF  at  High  Pressure 
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an  optical  cavity  composed  of  100%  and  70%  reflectors,  separated  by  1.  5 M 
was  used  over  the  entire  parametric  range. 

Variation  of  output  energy  with  partial  pressure  is  summarized  in 
Figures  71  and  72  for  total  pressures  of  1.7  and  2.3  atm. , respectively. 

The  variatior  In  output  energy  was  20-30%  due  to  shot-to-shot  repeatability. 

The  abscissa  in  each  p'  •'t  is  In  terms  of  capacitor  bank  charge  voltage  in 
order  to  provide  a consistent  means  of  comparison.  The  relationship  of 
this  parameter  to  actual  discharge  conditi&xis  is  presented  in  a subsequent 
graph.  For  both  sets  of  data,  the  Kr  partial  pressure  was  100  Torr,  considered 
to  be  optimum.  Variation  of  the  Kr  content  from  50-200  Torr  for  all  F^ 
partial  pressures  and  the  tota.  pressures  indicated,  resulted  in  only  small 
variations  of  the  'output  energy  on  the  order  of  the  experimental  variability. 

The  data  of  Figures  71  and  72  indicate  a general  trend  to  higher  output 
energies  with  increasing  charge  voltage  for  all  F^  partial  pressures. 

Sin^jular  inconsistencies,  for  example  the  datum  at  40  KV  for  1.5  Torr  F^ 
in  Figure  71,  were  not  representative  of  average  values  and  should  be 
ignored.  It  was  found  that  the  consistently  highest  outputs  were  obtained  at 
2 Torr  F^  (indicated  by  the  solid  lines)  and  a total  pressure  of  1.  7 atm.  As 
indicated  earlier,  laser  output  waveforms  at  higher  pressures  sMwed 
erratic,  spiking  behavior  therefore  experimentation  in  this  regime  was  not 
extensive. 

To  indicate  general  discharge  behavior  for  the  1.  7 atm  case,  peak 
discharge  voltage  and  currerc  at  laser  output  termination  are  shown  as  a 
function  of  charge  voltage  in.  F jure  73,  The  lower,  solid  ’-nes  indicate 
discharge  voltage  related  to  the  left  ordinate,  and  the  dotted  lines  represent 
terminal  current  related  to  the  right  ordinate.  For  increasing  F^  partial 
pressure,  discharge  voltage  is  seen  to  increase  whereas  current  decreases, 
indicating  higher  impedance  as  expected. 
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Ar  - 2. 3 atm 
Kr-  lOOTorr 
Jg  - 10-12  A/cm^ 
A-K  ■ 6 cm 
Cl -133  uF 


Figure  72.  Relative  Output  as  Function  of  ?2  Pressure 
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Ar  - 1.7  atm 
Kr  - lOOTorr 
Jg- 10-12  A/cm^ 
A-K  “ 6 cm 
C,  - a33  MF 


Charge  Voltage  (kV) 


Figure  73.  Peak  Voltage  and  Terminal  Current  as  Function  of  F2  Partial  Pressure 
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Terminal  Current  (KA)-Dotted  Lines 


As  previously  suggested,  a general  criterion  for  premature  discharge 

collapse  during  the  electron  gun  pulse  was  early  termination  of  the  laser 

output.  It  was  indicated  that,  above  a certain  critical  field,  higher  voltages 

would  lead  to  earlier  termination.  This  effect  is  shown  graphically  in  Figure 

74.  The  nominal  pulse  length,  with  electron  beam  excitation  only  was  800  ns 

as  indicated  at  the  abscissa  zero  point.  For  all  F partial  pressures, 

pulse  length  is  not  seriously  degraded  up  to  a charge  voltage  of  30  KV.  Above 

this  value,  however,  shortening  occurs  which  becomes  serious  for  higher  F^ 

partial  pressures.  For  2.0  Torr  of  F this  corresponds  to  a peak  E/N  of 
-17  2 ^ 

2.2  X 10  ' V-cm  . Of  course,  this  is  the  peak  experienced  during  the 

discharge  and  does  not  necessarily  correspond  to  the  undetermined  arcing 
value  which  occurs  later  in  the  pulse. 

Optimum  Results.  Laser  energy  extraction  was  optimized  for  both  small 
and  large  anode>cathode  separations.  In  the  first  case,  spacings  of  3-5  cm 
were  considered  in  the  hope  of  establishing  the  beam  deposition  uniformity 
required  for  high  enhancements.  As  will  be  shown  below,  discharge 
stability  is  crucially  dependent  on  deposition  uniformity  in  the  direction  of 
propagation.  Unfortunately,  the  results  in  this  instance  were  limited  by 
constraints  related  to  electrode  geometry.  In  the  second  case,  an  anode- 
cathode  sparl  g of  10  cm  was  used  in  the  interest  of  increased  energy  output. 

In  order  to  maximize  laser  output,  it  was  considered  important  to 
maximize  electron  beam  current  density.  In  this  regard,  the  70%  transmissive 
ground  plane  screen,  adjacent  to  the  gun  pressure  foil  was  removed.  At 
large  anode-cathode  spacings,  this  had  no  effect  on  discharge  stability,  but 
did  enhance  output  by  roughly  60%,  slightly  greater  than  would  be  expected 
by  the  simple  removal  of  a 30%  obscuration.  For  small  anode-cathode 
spacings,  however,  discharge  stability  was  greatly  impaired  as  a result  of 
arcing  to  the  raised  foil  edges  coincident  with  the  support  webs  which  become 
prominent  under  gas  pressure  loading.  In  order  to  provide  a smooth  ground 
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Laser  Pulse  Width  (ns) 


plane  without  reintroduction  of  the  70%  transmissive  screen,  a second 
-3 

unloaded  5 x 10  cm  aluminum  foil  was  placed  directly  over  the  pressure 
foil.  Measurements  with  a Faraday  cup  indicated  that  current  densities 
immediately  beyond  the  second  foil  were  identical  to  those  without  it.  Results 
with  this  configuration,  however,  were  not  significantly  different  from 
those  in  the  previous  arrangement  with  the  ground  plane  screen.  Possibly 
increased  beam  scattering  from  the  second  foil  was  responsible  for 
reducing  the  laser  output  below  expected  levels.  Alternatively,  attempts 
at  maximizing  enhancement  were  made  with  short  pulse  operation  at  an 
anode-cathode  separation  of  10  cm.  In  this  instance,  although  discharge 
collapse  was  experienced,  harder  discharge  pumping  could  be  achieved  over 
a brief  initial  period.  Therefore,  optimum  results  were  obtained  for  large 
volume,  short  and  long  pulse  excitation  as  outlined  below. 

Table  1 summarizes  the  results  for  long  pulse  operation  with  central 

2 

beam  current  densities  on  the  order  of  15-18  A/cm  . As  suggested  by 
previous  remarks  concerning  deposition  uniformity,  maximum  specific 
energy  was  deposited  roughly  in  the  central  three  liters.  Therefore,  in 
order  to  obtain  an  indication  of  maximum  specific  extraction,  the  laser 
output  was  masked  off  allowing  for  transmission  of  only  the  central  three 
liters  of  gain  volume.  The  three  liter  volume  was  calculated  assuming  an 
effective  gain  length  of  75  cm  as  indicated  by  the  longitudinal  deposition 
profiles.  These  results  are  tabulated  in  the  first  column  in  Table  1.  Specific 
extraction  in  this  case  was  relatively  high  at  18-Z2  J/liter.  The  span 
indicates  extraction  with  beam  only  for  the  lower  value  and  beam  plus 
discharge  excitation  at  the  maximum.  Enhancement  was  therefore  1.2. 

For  the  full,  unmasked  10  liter  extraction  volume  where  the  current 
gradient  was  estimated  at  10-18  A/cm^,  specific  extraction  was  7.  5-9.  0 
J /liter. 
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TABLE  1.  LARGE  SCALE  EXTRACTION  RESULTS  - LONG  PULSE 


3 liters 

10  liters 

Beam  Current 

15—18  A/cm^ 

E 

< 

oo 

f 

o 

Output  Energy 

54-65 J 

75-90  J 

Specific  Output 

18-22  J/C 

7.5-9.0J/C 

Specific  Input 

167-200  J/C 

100-120  J/C 

Efficiency 

11% 

7.5% 

Theoretical 

15  j/e 

Enhancement- 1. 2 

(ri  ■ 9%) 

Conditions 

Beam  Voltage:  320  kV 

Pulse  Length:  900  ns 

Optics:  R ■ 50%,  100% 

Gas:  Ar/kr/F2-  1290flOO/2,  1.7  atm 
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Enhancement  in  the  long  pulse  case  was  clearly  limited  by  premature 
discharge  collapse.  Therefore,  in  an  effort  to  increase  enhancement, 
harder  discharge  pumping  over  shorter  time  was  attempted.  The  relevant 
waveforms  are  shovm  in  Figure  75  and  the  results  are  summarized  in  Table 
2.  For  the  nominal  400  ns  pulse  length,  enhancement  was  increased  to  1.  8 
with  a specific  energy  output  of  7.2  J/liter.  As  indicated  in  the  second 
column  of  Table  2,  theoretical  predictions,  based  on  the  extensive  code 
described  earlier,  correlate  quite  well. 

Discharge  Power  Loading  Analysis.  As  indicated  above,  electron  beam 

deposition  uniformity  in  the  direction  of  propagation  was  the  limiting  factor 

In  achieving  high  discharge  enhancements.  A summary  of  a comprehensive 

analysis  is  presented  below  which  attempts  to  explain  these  limitations  in 

(58) 

terms  of  detailed  system  parameters'  . In  essence,  maximum 
discharge  energy  loading  is  limited  by  a local  stability  criterion  and  the 
requirement  for  current  continuity  from  anode  to  cathode. 

The  stability  criterion  is  derived  assuming  that  electron  production  and 
loss  are  a result  of  electron  beam  and  metastable  ionization  and  F^ 
attachment,  respectively.  An  upper  bound  on  discharge  energy  loading  is 
developed  from  a limitation  on  discharge  E/N,  which  if  beyond  certain  critical 
values  leads  to  runaway  metastable  ionization  and  arcing.  The  criterion 
reflects  the  requirement  that  ionization  from  the  electron  beam  must  dominate 
that  from  the  discharge  to  maintain  a stable  discharge  for  times  on  the  order 
of  a microsecond.  Or,  in  more  definite  terms,  for  the  parametric  range  of 
interest. 


(58) 

W.  H.  Long,  Jr,  "Discharge  Stability  and  Scaling  of  the  KrF  Laser", 
Presented  at  30th  Gaseous  Electronics  Conference,  Palo  Alto,  Calif. , 
x8-21  October  77,  (To  be  submitted  to  J.  Appl.  Phys.  (1977)). 
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TABLE  Z LARGE  SCALE  EXTRACTION  RESULTS  - SHORT  PULSE 


Experimental  Theoretical 

E-Beam  Only  E{ln)-83J/C  83J/C 

EIOut)-4J/)e  5.0J/C 

n ■ 5%  6% 

E-Beam  and  Discharge  Elln)  - 103  J/?  167  J/C 

ElOut)  - 7. 2 J/C  10  J/je 

n - 7%  6% 

Enhancement  ■ 7. 2/4  ■ L 8 2 

Conditions  E-Beam:  J ■ 18  A/cm^,  V ■ 300  keV,  T ■ 400  ns 

Discharge:  E/n  • 6 x 10'^^  V-cm^ 

Optics:  R ■ 5(W>,  100% 

Gas:  Ar/kr/F2  ■ 1290/100/2,  1 7 atm 
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(36) 


F < Y 

d e 


where,  P,  and  P are  discharge  and  electron  beam  power  density, 
d e 

respectively  and  the  quantity  y depends  on  only  the  fluorine  partial  pressure 
(F^)  and  total  pressure  (P)  according  to 

Y = 0.23  (F^)^’ Np)°*^  . (37) 

where  F^  is  in  Torr  and  P in  atm. 


It  is  required  that  the  current  be  continuous  across  the  discharge.  That 
is. 


J , = e n V = Const.  (38) 

d e 

Within  the  range  ITd  < E/N  < 10  Td  and  10  ^ < n/N  < 10  can  be 

approximated  by 

V = b(n/N)^’^  (E/N)°*^  . (39) 

e 

The  current  density  is  therefore, 

J = enb  (n/N)°*^  (E/N)°’^  (40) 

d 

_3 

The  requirement  (38)  therefore  dictates  that  E/N  vary  as  (n/N)  ; and  since 

P,  is  directly  proportional  to  E/N, 
d 

Pj  = (n/N)"^  . (41) 

d 

The  principal  loss  of  electrons  is  by  F^  attachment.  If  metostable  ionization 

is  neglected,  the  electron  density  can  be  expressed  as. 
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n = P /p  F,W. 
e 2 1 


(42) 


where  P is  the  attachment  rate  and  W.  is  the  energy  required  to  form  an 
electron-ion  pair. 

Combining  (41)  and  (42), 

-3 

P cc  P (43) 

d e 

For  the  purposes  of  calculation,  it  is  assumed  that  the  electron  energy 
deposition  fall-off  is  linear  from  the  foil  to  some  distance,  d,  from  the  foil. 

If  r is  the  ratio  of  power  density  at  d to  that  at  the  foil,  then 

P (x)  = P (l  + (r-l)x/d) 
e eo 

Therefore,  from  (43), 

P^(x)  = c (1  + (r-1)  x/d)’^  (44) 

d 

By  invoking  the  stability  criterion  (36)  and  choosing  the  constant  c to 
satisfy  this  at  the  point  d,  which  guarantees  it  for  x < d,  the  average 
discharge  power  loading  can  be  calculated  by  integrating  over  0 S x s d. 
Additionally,  if  a discharge  circuit  is  chosen  which  approximated  that  used 
experimentally,  then  the  resulting  power  loading  must  be  divided  by  2,  due 
to  the  approximately  sinusoidal  waveform. 

Finally, 


P = P (l  + r)/2 

(45) 

e eo 

P = Y P^^^r^^^  (l+r)/2 
d eo 

(46) 
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The  enhancement,  c,  defined  by 


e 


P + P 
d e 

P 

e 


is  therefore  given  by. 


e 


eo 


(47) 


This  relationship  is  graphed  in  Figure  76  for  a beam  power  density  of  0.4 

3 2 

MW/cm  , corresponding  to  a current  density  of  15  A /cm  , and  a nominal 

mixture  of  F^/Kr/Ar  = 2/100/1267  Torr.  For  the  range  of  enhancements 

obtained  experimentally,  (1. 2-1.4),  the  indicated  value  of  r is  about  1/3 

and  this  is  well  within  the  values  of  beam  deposition  determined  experimentally 

as  outlined  in  Section  2.  Also,  by  (47),  enhancement  decreases  as  beam 

current  is  raised,  and  this  was  observed  experimentally.  In  this  case,  the 

theoretical  maximum  enhancement  is  only  1.9.  Of  course,  for  fixed  r, 

enhancement  can  be  increased  by  increasing  F^  partial  pressure  and  total 

pressure  or  by  reducing  electron  beam  power  density.  However,  these 

procedures  are  limited  due  to  kinetic  processes.  As  beam  current  density 

is  reduced,  total  power  drops  and  absorption  begins  to  become  important. 

.JU 

For  an  increase  in  F^,  quenching  of  KrF  and  absorption  by  F^  increases; 
and  an  increase  in  total  pressure  accentuates  beam  nonuniformity  and  leads 
to  the  formation  of  timers  which  provide  additional  energy  loss  channels. 
Determination  of  the  optimum  set  of  parameters  is  very  complex  and  cannot 
^e  obtained  simply  from  the  above  formulation. 


4.  2 Typical  Theoretical  Correlation.  The  details  of  an  extensive 
theoretical  kinetic  model  for  description  of  the  physics  of  an  electrically 
excited  laser,  as  well  as  its  numerical  implementation  by  a comprehensive 
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and  general  computer  code,  have  been  described  in  Section  3.  In  Section 
3.4,  results  of  a typical  analysis  were  presented  in  order  to  illustrate  the 
flexibility  and  generality  of  this  code;  in  addition,  the  figures  presented  in 
that  Section  have  been  selected  to  summarize  the  physics  of  the  molecular 
kinetic  reaction  scheme  currently  assumed  for  the  KrF  laser  system,  as  well 
as  to  correspond  to  parameters  for  which  a comparison  with  experimental 
data  from  the  scale-up  device  can  be  made.  Correlations  of  the  theory  with 
experimental  data  have  been  made  for  a variety  of  conditions,  and  reasonably 
good  agreement  has  been  obtained. 

The  nominal  operating  parameters  for  the  calculations  presented  in 

Section  3.4  correspond  to  conditions  typical  of  the  scale-up  device  described 

earlier.  This  device  consists  of  an  electron-beam-stabilized,  low  inductance, 

capacitive  discharge,  defined  by  a 20  x 100  cm^  anode,  separated  by  (its 

maximum  value)  of  10  cm  from  a grounded  cathode  screen  through  which  a 

(~  320  keV)  e -beam  is  passed.  The  maximum  g.-m  volume  is  nominally 

20  liter,  with  optical  extraction  along  the  major  axis.  However,  for  the 

present  experiment,  the  optical  resonator  used  a 6 cri  diameter  (R  = 70%, 

not  optimized)  extraction  area  located  just  adjacent  to  the  screen  cathode 

with  L (cavity  spacing)  = 130  cm,  and  y (the  intracavity  loss/pass)  = 0. 
c 

The  e -beam-excited  gain  length  was  estimated  (by  cinemoid)  to  be  L ~ 75  cm, 

S 

giving  an  extraction  volume  ~ 2 liter.  The  laser  plenum  can  be  operated  up 
to  10  atm,  although  the  typical  conditions  for  the  present  results  consisted  of 
Ar/Kr/F^  = 1190/100/1.  5 at  a total  pressure  of  1.7  atm.  Electrical  pumping 
was  provided  by  direct  e -beam  exc  tion,  or  enhance.*  by  an  electrical 
discharge  with  energy  supplied  by  a 1.  33  pF  capacitor  bank,  initially  charged 
to  40  kV.  Measured  circuit  inductance  was  200  nH.  Electron  beam  current  density 
was  ~ 10-12  A/cm^  for  a pulse  duration  of  < 1.0  us.  Available  circuit 
diagnostics  included  discharge  voltage,  measured  directly  at  the  anode,  and 
total  discharge  current,  measured  at  the  capacitor  bank  ground  connection. 

Optical  output  was  diagnosed  with  a photodiode  for  temporal  behavior,  and  a 
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calorimeter  for  measurement  of  total  output  energy.  Predictions  of  the 

analytical  model  have  been  compared  with  experimental  results  from  this 

device.  Figure  77a  shows  an  oscilloscope  trace  of  the  e -beam  voltage 

and  Figure  77b  displays  the  temporal  output  power  obtained  with  e -beam 

excitation  only.  Figure.  78a  shows  traces  for  the  discharge  voltage  and 

current,  and  Figure  78b  the  laser  output  enhanced  with  discharge  pumping. 

These  data  are  to  be  compared  with  theoretical  predictions  shown  in  Figure 

3/2 

79.  The  e -beam  current  density  J,  (t),  which  approximately  follows  a V 

^ . 2 

law,  had  a maximum  value  of  ~ 12  A/cm  , and  is  shown  in  Figure  79. 

Experimentally,  the  maximum  extracted  output  was  ~8  J with  e -beam 
excitation  only,  and  14  J with  an  electric  discharge,  corresponding  to 
specific  outputs  of  4.0  J/liter  and  7.0  J/liter,  respectively.  Beam  uniformity 
was  estimated  to  be  approximately  constant,  with  less  than  15%  reduction  in 
intensity  at  the  edges  of  the  optical  aperture.  These  values  are  to  be  compared 
with  predictea  values  (integration  of  power  in  Figure  79)  ox  4.7  and  7.0  J/liter, 
at  5.6%  and  5.5%  efficiency,  respectively. 

Note  the  steep  leading  edge  and  the  sharp  decay  at  early  times  on  the 

optical  extraction.  The  origin  of  this  effect  is  the  theoretically  predicted 

strong  overshoot  of  the  net  gain  '.oefficient  (about  a factor  of  5)  over  the 

threshold  loss,  resulting  in  superfluorescent  extraction  for  about  100  ns 

(cavity  photon  lifetime  is  t^  ~ 25  ns),  after  which  the  gain  relaxes  to 

threshold.  Later  in  the  pulse,  for  e -beam  excitation  only,  the  output 

approximately  follows  the  temporal  e -beam  current  density,  sine®  the 

dominant  mechanism  is  the  formation  of  KrF  by  F ion  recombination. 

(The  F population  density  follows  the  secondary  electron  density  n , whose 

e 

source  term  is  e -beam  ionization,  J,  (t)a  N/e).  For  the  case  of  electric 

b eff 

discharge,  the  rate  of  change  of  sustainer  voltage  and  current  is  limited  by 
the  inductance,  and  the  maximum  voltage  ("-  12.  5 kV)  is  limited  by  the  low 
circuit  impedance  (^  0.2  Cl).  Thus,  the  effects  of  enhanced  pumping  do  not 
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J ' 200  ns 


2.5  kV/div 
18  kA/div 


7 J/liter 


E Figure  78.  Discharge  Characteristics  and  Optical  Output  Power 

p 
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Figure  79.  Thooretically  Predicted  Optical  Output  Power  Density  With/Without  Discharge 


occur  until  much  later  in  the  pulse,  ifter  ~ 500  ns.  Note  the  appearance  in 
Figure  78a  of  a sharp  initial  voltage  spike  across  the  discharge,  caused  by 
the  initially  high  impedance.  Note  also  the  increase  in  discharge  voltage 
at  the  end  of  the  pulse,  after  the  e -beam  current  falls  off.  This  is  also 
reflected  in  Figure  78b  and  79  in  the  small  bump  on  the  final  edge  of  the 
optical  extraction  for  the  discharge  case  (somewhat  overpredicted  by  the 
analysis).  This  perturbation  occurs  because  the  discharge  impedance  rises 
faster  than  the  current  decays. 

The  good  agreement  which  has  been  obtained  in  comparisons  of 
predictions  of  the  present  model  with  data  from  this  and  other  KrF  laser 
devices  suggests  that  the  KrF  kinetics  are  reasonably  well  understood. 
Results  of  the  present  analysis  indicate  that  much  higher  optical  extraction 
should  be  possible  with  the  nominal  operating  parameters  of  the  present 
device.  The  extrapolation  to  higher  power  will  be  discussed  in  the  following 
section. 

4.  3 Scaling  Considerations.  Predictions  of  the  present  laser  kinetic 
model  indicate  that  high  optical  extraction  (>  35  J/liter)  at  ~ 6.5%  efficiency 
should  be  possible  with  the  nominal  operating  parameters  of  the  present 
device.  To  achieve  this  extraction,  enhancement  by  means  of  an  electric 
discharge  is  required.  At  present,  the  optimistic  performance  calculations 
to  be  discussed  here  have  not  been  realized  experimentally^  and  the  reason 
for  this  can  be  explained  in  terms  of  discharge  nonuniformity  and  stability. 

If  uniform  e -beam  energy  deposition  and  uniform  electrical  excitation  over 
large  volumes  can  be  achieved,  the  realization  of  efficient  high  energy  KrF 
laser  performance  should  be  possible. 

Consider  the  following  experimental  parameters  (optimistic,  but  not 
necessarily  optimized)  which  shoulc  be  attainable  from  the  present  scale-up 
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device.  Assume  an  Ar/Kr/F^  = 89.8/10.0/0.2  mixture  at  3.  0 atm,  pumped 

for  1.  0 (i.  s with  an  external  RLC  circuit  with  L = 200  nH  and  C = 4.  0 pF 

2 

(initially  charged  to  100  kV)  driving  a discharge  of  area  2000  cm  and  A-K 

separation  of  10  cm,  sustained  by  an  external  e -beam  of  (constant)  current 

density  15  A/cm^.  Let  the  optical  resonator  be  formed  by  (100,  30)%  reflecting 

mirrors  with  a gain  length  L = 100  cm.  Results  of  the  analysis  are  shown  in 

S 

Figure  80,  where  optical  output  power  density  and  electrical  input  power 
densities  from  the  e-beam  and  discharge  are  shown  plotted  as  a function  of 
time.  For  a 1.0  ps  pulse,  integration  of  the  curves  gives  an  output  energy 
extraction  of  > 35  J/liter  at  ~ 6.5%  efficiency. 

4.4  Conclusions.  Results  of  the  present  investigation  have  indicated 
that  high  optical  energy  extraction  at  high  efficiencies  is  possible  from  the 
KrF  system.  As  a result  of  the  present  program,  an  extremely  flexible 
and  versatile  experimental  device,  as  well  as  a comprehensive  theoretical 
analysis  with  a broad  range  of  applicability  were  developed  for  the  investigation 
of  electrically  excited  laser  systems.  Although  certain  experimental  limitations 
prevented  the  realization  of  initial  program  goals  for  the  KrF  laser  system, 
significant  results  were  achieved  and  understanding  of  fundamental  physical 
mechanisms  was  obtained. 
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